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We use both Mutual Composite Fermion (MCF) and Composite Boson (CB) approach to study 
balanced and im-balanced Bi-Layer Quantum Hall systems (BLQH) and make critical comparisons 
between the two approaches. We find the CB approach is superior to the MCF approach in studying 
ground states with different kinds of broken symmetries. In the phase representation of the CB 
theory, we first study the Excitonic superfluid state (ESF). The theory puts spin and charge degree 
freedoms in the same footing, explicitly bring out the spin-charge connection and classify all the 
possible excitations in a systematic way. Then in the dual density representation of the CB theory, 
we study possible intermediate phases as the distance increases. We propose there are two critical 
distances d c i < al C 2 and three phases as the distance increases. When < d < d c i, the system is in 
the ESF state which breaks the internal U(l) symmetry, when d c i < d < dci, the system is in an 
Pseudo-spin density wave ( PSDW ) state which breaks the translational symmetry, there is a first 
order transition at d c i driven by the collapsing of magneto-roton minimum at a finite wavevector 
in the pseudo-spin channel. When d C 2 < d < oo, the system becomes two weakly coupled v = 1/2 
Composite Fermion Fermi Liquid ( FL) state. There is also a first order transition at d — d C 2- We 
construct a quantum Ginzburg Landau action to describe the transition from ESF to PSDW which 
break the two completely different symmetries. By using the QGL action, we explicitly show that 
the PSDW takes a square lattice and analyze in detail the properties of the PSDW at zero and finite 
temperature. We also suggest that the correlated hopping of vacancies in the active and passive 
layers in the PSDW state leads to very large and temperature dependent drag consistent with the 
experimental data. Then we study the effects of imbalance on both ESF and PSDW. In the ESF 
side, the system supports continuously changing fractional charges as the imbalance changes. In the 
PSDW side, there are two quantum phase transitions from the commensurate excitonic solid to an 
in-commensurate excitonic solid and then to the excitonic superfluid state. We also comment on the 
effects of disorders and compare our results with the previous work. The very rich and interesting 
phases and phase transitions in the pseudo-spin channel in the BLQH is quite similar to those in 4 He 
system with the distance playing the role of the pressure. A BLQH system in a periodic potential is 
also discussed. The Quantum Hall state to Wigner crystal transition in single layer Quantum Hall 
system is studied. 



I. INTRODUCTION 

(1 ) Experimental Observations 

Extensive attention has been lavished on Fractional 
Quantum Hall Effects (FQHE) in multicomponent sys- 
tems since the pioneering work by Halperin 1 . These com- 
ponents could be the spins of electrons when the Zeeman 
coupling is very small or layer indices in multi-layered 
system. In particular, spin-polarized Bilayer Quantum 
Hall systems at total filling factor vt — 1 have been un- 
der enormous experimental and theoretical investigations 
over the last decade^. When the interlayer separation 
d is sufficiently large, the bilayer system decouples into 
two separate compressible v = 1/2 layers. Earlier ex- 
periments exhibited a strong suppression of the tunnel- 
ing current at low biases^. However, when d is smaller 
than a critical distance d c , even in the absence of inter- 
layer tunneling, the system undergoes a quantum phase 
transition into a novel spontaneous interlayer coherent 
incompressible phased. At low temperature, with ex- 
tremely small interlayer tunneling amplitude, Spielman 
et al discovered a very pronounced narrow zero bias peak 
in this interlayer coherent incompressible stated. M. Kel- 



logg et al also observed quantized Hall drag resistance at 
h/e 2 —. In recent counterflow experiments, it was found 
that both linear longitudinal and Hall resistances take 
activated forms and vanish only in the low temperature 
limits. 

(2) Theoretical achievements 

Starting from Halperin's 111 wavefunction Eqn IAll 
which describes a bi-layer system with Ni ( N2 ) electrons 
in the top (bottom) layer ( the total number of electrons 
is N = Ni + N2 ), using various methods, several authors 
discovered a Neutral Gapless Mode (NGM) with linear 
dispersion relation uj ~ vk and that there is a finite tem- 
perature Kosterlitz-Thouless (KT) phase transition as- 
sociated with this NGM&SJii. By treating the two layer 
indices as two pseudo-spin indices, Girvin, Macdonald 
and collaborators mapped the bilayer system into a Easy 
Plane Quantum Ferromagnet (EPQFM) ^ 11 ' 12 ( which is 
equivalent to the Excitonic Superfluid ). They estab- 
lished the mapping by projecting the Hamiltonian of the 
BLQH onto the Lowest Landau Level (LLL) and then 
using subsequent Hartree-Fock (HF) approximation and 
gradient expansion ( called LLL+HF in the following ). 
In the picture of EPQFM^, the canonical ensemble with 
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definite S z = N\ — N2 = is replaced by Grand canon- 
ical ensemble with fluctuating S z . The relative fluctu- 
ation of S z is at the order of 1/y/N — ► as N — ► 00. 
By drawing the analogy with superconductivity where 
canonical ensemble with definite number of Cooper pairs 
is shown to be equivalent to BCS wavefunction which 
is a grand canonical ensemble with indefinite number of 
Cooper pairs, the authors in£ i n i 12 argued this trial wave- 
function is a good approximation to the exact ground 
state. The low energy excitations above the ground state 
is given by an effective 2 + 1 dimensional XY model. 
There are 4 flavors of topological defects called " merons 
" which carry fractional charges ±1/2 and also have ± 
vorticities. They have logarithmic divergent self energies 
and are bound into pairs at low temperature. The lowest 
energy excitations carry charge ±e which are a meron 
pair with opposite vorticity, but the same charge. There 
is a finite temperature phase transition at Tkt where 
bound states of the 4 flavors of merons are broken into 
free merons. The large longitudinal resistivity (~ lfcfi) 
observed ini at very low temperature indicated that these 
meron pairs are highly mobil. In the presence of small 
tunneling, thsy^L found that when the applied in-plane 
magnetic field is larger than a critical field Sm, there is 
a phase transition from a commensurate state to an in- 
commensurate state ( C-IC) with broken translational 
symmetries. When B > B^, there is a finite temperature 
KT transition which restores the translation symmetry 
by means of dislocations in the domain wall structure in 
the incommensurate phase. Starting from the EPQFM 
approach, several groups investigated I — V curves in the 
presence of small tunneling—. In addition to the work 
mentioned above, there are also many other works done 
on BLQH. For example, several authors applied differ- 
ent versions of composite fermion Chern-Simon theory 
to study BLQH systems in 16 i 17 i 19 ,. 

(3) Discrepancies between theory and the experi- 
ments. 

Despite the intensive theoretical research 
in^&ii i 10 ! 11 ! 12 , there are still many serious discrep- 
ancies between theory and the experiments. According 
to the present theories, in the Excitonic Superfluid ( 
ESF ) state, there should be the interlayer tunneling 
Josephson effect, finite temperature KT transition and 
vanishing linear longitudinal and Hall resistances in the 
counterfiow channel below the KT transition tempera- 
ture. Unfortunately, all these characteristics have never 
been observed in the experiments. Although it appears 
certain that the dramatic conductance peak observed at 
d < d c is due to the collective NGM in the ESF state, 
no theory can explain the magnitude of zero bias peak 
conductance which, though enormously enhanced, does 
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not exceed 10 its width and its dependence on 

proximity to the ESF state boundary ( see, however, 
the most recent work-i^ ). All the previous calculations 
predicted that the in-plane field will split the zero bias 
peak into two side peaks, but the experiments showed 
that although there are two tiny shoulders appearing, 



the central peak stay. The excess dissipation observed 
in the counterfiow experiments^ appears to vanish only 
as T — > limit. The origin of the excess dissipation 
remains an important unresolved question. But it was 
argued ini that the origin may come from disorders 
induced mobil vortices. 

(4) Open problems in SLQH and BLQH 

One of the remaining outstanding problems in single 
layer quantum Hall ( SLQH ) system is to understand 
the quantum phase transitions from QH to QH or QH to 
insulating state as tuning the magnetic field at a fixed 
electron density or vice versa. Similarly, one of the re- 
maining outstanding problems in BLQH is to understand 
novel phases and quantum phase transitions as changing 
the distance between the two layers. When the distance 
is smaller than d c \, the system is in the ESF state, while 
it is sufficiently large, the system becomes two weakly 
coupled v — 1/2 Composite Fermion Fermi Liquid (FL) 
state. There could be a direct first order transition 
between the two states as indicated in some numerical 
calculations 1 --. However, the experimental observations 
that both zero voltage tunneling peak* and the Hall drag 
resistivity^ develop very gradually when d ~ d c \ suggest 
the transition at d = d c \ is a 2nd order phase transi- 
tion. Although there are very little dissipations in both 
the ESF and FL, the experiment 2 ^ discovered strong en- 
hancement of drag and dissipations in a large interme- 
diate distance range. These experimental observations 
suggest that there must be intermediate phases between 
the two phases. Unfortunately, so far, the nature of the 
intermediate phase, especially the quantum phase tran- 
sition between the ESF and the intermediate phase was 
not systematically investigated. The outstanding prob- 
lems in both SLQH and BLQH will be addressed in this 
paper. 

(5) Advances in Helium 4 system 

A superfluid is a fluid that can flow through the tiniest 
channels or cracks without viscosity. The phenomenon of 
superfluidity was first observed in the two isotopes of He- 
lium: i He and 3 He which become superfluids below the 
transition temperatures T c = 2.18K and T c — 2AmK 
respectivel y 21 ! 22 . Recently, a PSU group lead by Chan 
found some signature of a possible supersolid state in 
A He at low temperature T < 200 mK and high pres- 
sure p > p c ~ 25 bar. The authors i n 25 i 26 ' 27 argued 
that there may be vacancies even at T = whose con- 
densation leads to the supersolid 4 He. In2£, the author- 
constructed a Ginsburg Landau (GL ) theory to map out 
the He phase diagram, analyze carefully the conditions 
for the existence of the supersolid (SS) and study all the 
phases and phase transitions in a unified framework. I 
introduced a single parameter g which is the coupling 
between the normal solid (NS)component and superfluid 
(SF) component in the GL theory. If g = g v < 0, I ex- 
plicitly showed that there are two scenarios (1) If \g v \ 
is sufficiently small, then the normal solid is a commen- 
surate solid ( C-NS ). The C-NS still does not have a 
particle-hole symmetry, it is a vacancy-like NS where the 
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excitation energy of a vacancy is lower than that of an in- 
terstitial, therefore named NS-v. Then for the first time, 
we construct a quantum Ginsburg-Landau theory (QGL) 
to study the SF to NS-v transition and explicitly derive 
the Fcymann relation from the QGL (2) If \g v \ is suf- 
ficiently large, then a vacancy induced SS ( SS-v ) ex- 
ists at sufficient low temperature. The critical tempera- 
ture Tss-v becomes an effective measure of the coupling 
strength g v . The analogy between ESF in the pseudospin 
channel in the BLQH and the superfluid in A He was ex- 
plored by many previous work^ i 8 ' 9 ' 10 ? 11 ! 12 . In this paper, 
we push this analogy to deeper and broader contexts. We 
find that the very rich and interesting phases and phase 
transitions in the pseudo-spin channel in the BLQH is 
quite similar to those in A He system with the distance 
playing the role of the pressure. 

(6) The analogy between BLQH and Helium 4 to be 
explored in this paper 

In this paper, we use both Mutual Composite 
Fermion ( MCF^iM^ and Composite Boson ( CB) 
approaches 37,38 to study balanced and im-balanced 
BLQH systems. We identify many problems with MCF 
approach. Then we develop a simple and effective CB 
theory which naturally avoids all the problems suffered in 
the MCF approach. The CB theory not only can describe 
many crucial properties of the ESF elegantly, it can also 
be applied to study the possible novel intermediate phase 
with broken translational symmetry and the quantum 
phase transitions between these different ground states. 
Therefore the CB approach is superior to the MCF ap- 
proach in describing the BLQH where ground states with 
different kinds of broken symmetries may happen as the 
distance changes. By using this CB theory in its phase 
representation, we first study the excitonic superfluid ( 
ESF ) state. The theory naturally puts spin and charge 
degree freedoms in the same footing, explicitly bring out 
the spin-charge connection and classify all the possible 
excitations in a systematic way. Then by using the CB 
theory in its dual density representation and inspired by 
the insights achieved from the QGL theory developed in^ 8 . 
to describe the superfluid to solid and supersolid transi- 
tion in 4 He, we explore many similarities and also some 
difference between the BLQH and the A He system. Our 
understanding on phases and phase transitions in A He 
system can shed considerably lights on the phase and 
phase transitions in BLQH with the pressure playing the 
role of the distance. For balanced BLQH, there are two 
critical distances d c \ < d c2 . When < d < d&, the sys- 
tem is in the ESF state, when d c \ < d < d C 2, the system 
is in the Pseudo-spin density wave (PSDW) state, there 
is a first order transition at d c \ driven by magneto-roton 
minimum collapsing at a finite wavevector in the pseudo- 
spin channel. When d C 2 < d < oo, the PSDW melts into 
two weakly coupled v = 1/2 Composite Fermion Fermi 
Liquid (CFFL) state. There is a also first order transition 
at d = d C 2- However, disorders could smear the two first 
order transitions into two second order transitions. The 
transition from the ESF to the PSDW is unusual because 



the two states break two completely different symmetries: 
the global internal U(l) symmetry and the translational 
symmetry respectively. We construct an effective theory 
in the dual density representation to describe this novel 
quantum phase transition. We find the transition is very 
similar to the superfluid to normal solid transition in He 
system with the distance playing the role of the pressure. 

(7) Experimental setup 

Consider a bi-layer system with N± ( N2 ) electrons in 
left ( right ) layer and with interlayer distance d in the 
presence of magnetic field B = V x A ( Fig.l): 

H = H + Hi n t 

H = d x& a (x) — - c a {x) 

J 2m 

H int = ^ J d 2 xd 2 x'8p a (x)V a(i {x - x' )Spp(x') (1) 

where electrons have bare mass m and carry charge — e, 
c a ,a = 1,2 are electron operators in top and bottom 
layers, 8p a (x) — c^x^c^x) — n a ,a = 1,2 are normal 
ordered electron densities on each layer. The intralayer 
interactions are V\\ = V22 = e 2 /er, while interlayer in- 
teraction is V12 = V21 — e 2 /e\/r 2 + d 2 where e is the 
dielectric constant. For imbalanced bi-layers, ri\ 7^ n 2 , 
but the background positive charges are still the same 
in the two layers, the chemical potential term is already 
included in H int in EqnfT] There are two limits: (1) 
Weak tunneling limit A S as <C — • (2) Strong tunnel- 
ing limit Asas 3> — where the bilayer system becomes 
essentially the same as a single layer. In this paper, we 
only consider the non-trivial weak tunneling limit at to- 
tal filling factor vt = v\ + v 2 = 1 . For simplicity, we set 
&SAS = in the EqnQ] 

(8) The organization of the paper 

The rest of the paper is organized as follows. In section 
II, we use a MCF theory to study the BLQH. We achieve 
some limited success, but also run into many troublesome 
problems. In section III, we use a CB theory which puts 
spin and charge sector on the same footing to study the 
BLQH. We demonstrate why this CB theory naturally 
avoid all the problems suffered in the MCF approach. 
We also compare the CB approach with the EPQFM 
approach and point out advantages and limitations of 
both approaches. In section IV which is the key section 
of this paper, we push the CB theory further to ana- 
lyze carefully the instability in the pseudo-spin channel 
which lead to the PSDW state and study its properties 
at both zero and finite temperatures. We also point out 
the important physical consequences of possible quantum 
fluctuation generated vacancies in the PSDW state. We 
also compare our results with previous ones on the in- 
termediate phase. We study the effects of imbalance on 
both ESF and PSDW in the corresponding sections III 
and IV. In the final section, we summarized the main 
results of the paper and concluded that CB theory is su- 
perior to MCF approach in BLQH systems with ground 
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FIG. 1: A fabricated bilayer quantum Hall system. The 
quantum well width w, the well center to well center distance 
d and the magnetic length I ( which is the average distance 
between electrons in each well) are all of the same order ~ 
100A. Asas is the tunneling between the two wells. 



states of different broken symmetries. In Appendix A, 
we evaluate the meron fractional charges from its trial 
wavefunction in imbalanced case. In Appendix B, we in- 
vestigate one of the outstanding problems in SLQH: the 
Quantum Hall state to Wigner crystal state transition in 
SLQH. In appendix C, we study phases and phase tran- 
sitions in a BLQH in a periodic potential and two other 
closely related interesting models. 



II. MUTUAL COMPOSITE FERMION 
APPROACH: LIMITED SUCCESS AND FAILURE 

In parallel to advances in bi-layer QH systems, much 
progress has been made on novel physics involving quasi- 
particles and vortices in high temperature superconduc- 
tors. Anderson employed a single- valued singular gauge 
transformation to study the quasi-particle energy spec- 
trum in the vortex lattice stated. By employing the 
Anderson transformation, the author studied the quasi- 
particle transport in random vortex array in the mixed 
state 30 . The author also extended the Anderson singu- 
lar gauge transformation for static vortices to a mutual 
singular gauge transformation for quantum fluctuation 
generated dynamic vortices 3 -^. By using this dynamic 
gauge transformation, the author investigated the zero 
temperature quantum phase transition from d-wave su- 
perconductor to underdoped side by assuming the tran- 
sition is driven by the condensations of quantum fluctu- 
ation generated vortices 3 -^. In this section, by employing 



essentially the same singular gauge transformation used 
to study the interactions between quasi-particles and vor- 
tices in high temperature superconductors, we revisit the 
bi-layer QH systems. 

(1 ) Singular Gauge Transformation: 

Performing a single- valued singular gauge transforma- 
tion (SGT) 3 ^: 



U 



i(p J d 2 x J d 2 x' p\ (x)arg(x— x')p2 (%') 



i 



(2) 



we can transform the above Hamiltonian into: 



(3) 

where the transformed fermion is given by: 

V>i(£) = U Cl {x)U~ l = Cl {x)e^ d2x ' arg{s - 3 '^ 2 ^ 
i(j 2 (£) = Uc 2 {x)U- 1 = c 2 {x)e l ^ cpx ' arg{s '-^ pl ^{A) 

and the two mutual Chern-Simon (CS) gauge fields a a 
satisfies: V-a Q = 0,Vxo a = 2np a {x) ( See Fig. 2a). Ob- 
viously, the interaction term is unaffected by the singular- 
gauge transformation. Note that it is arg(x — x') ap- 
pearing in %j}\{x), while arg(x' — x) appearing in ip2{x) 
in EqnJU This subtle difference is crucial to prove all 
the commutation relations are kept intact by the single- 
valued singular gauge transformation EqnJ5] Note also 
that arg(x' — x) works equally well in Eqn[21 

It is easy to check that in single layer system where 
P\ = 92 = P, EqnJ5] reduces to the conventional singular 
gauge transformation employed in 35 : 



a i2 f d 2 x' arg{x — x)p a (x) 



Ca(x) 



(5) 



where p a (x) = c\(x)c a {x) is the electron density in layer 
a = 1, 2. It puts two flux quanta in the opposite direc- 
tion to the external magnetic field at the position of each 
electron ( Fig. 1 b). On the average, a CF feels a reduced 
effective field which is the external magnetic field minus 
the attached flux quanta. 

For i/ T = 1 bi-layer system, Eqn0]puts one flux quan- 
tum in one layer in the opposite direction to the external 
magnetic field at the position directly above or below 
each electron in the other layer ( Fig. 2 a). 

On the average, a Mutual Composite Fermion (MCF) 
in each layer feels a reduced effective field which is the 
external magnetic field minus the inserted flux quanta 
in this layer. In single layer system, it is essential to 
attach even number of flux quanta to keep Fermi statis- 
tics intact. The two attached flux quanta are moving 
together with the associated electron. However, in bi- 
layer system, inserting one flux quantum does keep Fermi 
statistics intact. The inserted one flux quantum in one 
layer is moving together with its associated electron in 
the other layer. If choosing 0=1/2, the transformation is 
not single- valued, the statistics is changed from fermion 



5 



B 




d<d c 

(a) 



d » d c 
(b) 



FIG. 2: Contrast the flux attachment in Eqn[2](a) with that 
in Eqn[S](b). In (a), there is one flux quantum in layer 1 when 
there is an electron directly downstairs in layer 2 and there is 
one flux quantum in layer 2 when there is an electron directly 
upstairs in layer 1. One can compare Fig. la with Fig. 2. 
For simplicity, we only show three electrons in the top layer 
and two electrons in the bottom layer. 



layer 2 will create a charge e 2 = — = — v 2 at x ( note 
that el+e 2 = —1 ). The product of operators ipl(x)?p2(y) 
in layer 2 will create a pair of fractional charges ±e 2 at x 
and y with the charge gap hiu* 2 . Only when the two lay- 
ers are identical iVi = N 2 , e\ = e 2 = —1/2 which carries 
fractional charge of even denominators. In general, the 
imbalanced bilayer system supports continuously chang- 
ing total fractional charges in the thermodynamic limit. 

The above arguments give the correct total fractional 
charges ±e{,±e 2 . However, it can not determine the 
relative charge distributions between the two layers. At 
mean field level, the energies of all the possible relative 
charge differences are degenerate. The lowest energy con- 
figuration can only be determined by fluctuations. The 
above arguments are at most intuitive. The much more 
rigorous and elegant topological arguments can only be 
given in the composite boson theory of the next section. 

(4) Fluctuations: 

When considering fluctuations around the MCF mean 
field theory, it is convenient to go to Lagrangian 36 : 



C = ipl{d T - iao)ip a + 



( — iV + A a - a a f 
2m% 



■Ipa 



to boson, this choice will be pursued in the next section 
on Composite boson approach. 
(2) Mean field theory: 

In the following, we put H = c = e = e = l. At total 
filling factor vt = 1, V x A = 27m where n = n\ + n 2 
is the total average electron density. By absorbing the 
average values of C-S gauge fields Vx < d a >= 2mi & 
into the external gauge potential A* a = A— < a a >, we 
have: 



H 



d 2 xip\,(x) 



HV + A* a (x) -6d a (x)f 
2m 



ip a {x) (6) 



where V x A* a = 2itn a and V x Sa a = 2-K5pa(x) are the 
deviations from the corresponding average density ( In 
the following, we will simply use a a to stand for these 
deviations). 

When d < d c , the strong inter-layer interactions renor- 
malize the bare mass into two effective masses m*— . 
MCF in each layer feel effective magnetic field -B* = 
V x A* a — 2iTn a , therefore fill exactly one MCF Lan- 
dau level. The energy gaps are simply the cyclotron gaps 
of the MCF Landau levels u>* = ^1. 

COL m * 

(3) Fractional charges: 

Let's look at the charge of quasi-particles created by 
MCF field operators ipa{x). Intuitively, inserting i/j\(x) 
on layer 1 not only inserts an electron in layer 1 at posi- 
tion x, but also pushes v 2 = ^ electrons in layer 2 into 
its boundary, therefore induces a local charge deficit at x 
in layer 2 which carries charge V2, the total charge at x is 
e\ = — 1 + v 2 = — v\. If x and y are two points far apart, 
then the product of operators ^\{x)^)i{y) in layer 1 will 
create a pair of fractional charges ±ej at positions x and 
y with the charge gap hw* cl . Similarly, inserting tp\ (x) in 



+ ia^na + — (a 1 ) a /3aQa'j: 



-in 



(a i t a i t +aiai + 2e-^a i t ai) 



(7) 



where the constraints have been used to rewrite the 
Coulomb interactions and af is the transverse spa- 
tial component of gauge field in Coulomb gauge V • 
a a = 35 ' 76 i 77 ( In Lorenz invariant gauge the sec- 
ond to the last term becomes the mutual C-S term 31 
j^{cFi)ap^tiu\o^d v a^ where a, (3 refer to layer indices, 
while fi, A refer to space-time indices ). The 
Hamiltonian has local x UiX)i gauge symmetry 

which corresponds to the invariance under ijj a (x) — > 



Integrating out MCF -0i , tfj 2 to one-loop and carefully 
expanding the interlayer Coulomb interaction to the nec- 
essary order in the long-wavelength limit leads to: 



2tt 



4tt 



(at 



+ 



^ 2 (a + ) 2 + i( e+ c 2 + ( X+ -^) ( z 2 )(a+) 2 

^ ? 2 («o ) 2 + -M" 2 + (X+ + ^)9 2 )K") 2 



27T' 



+ —q 2 a+a + -{e-Uj 2 + X-qW t a t + 



2 ' u u 2 

where • • • are higher gradient terms and af is the trans- 
verse component of gauge field in Coulomb gauge V-a a = 

°i = a l ± « 2 and e ± = 5( £ i ± e 2 ),X± = 
stands for the total ( rela- 



( ^ ) 



n35,76,77 „± 

^™ — ; <v 

\{Xi ± Ac- 
tive ) density fluctuation. a~ is the NGM identified 

previously&2i±&. The dielectric constants e a = and 

the susceptibilities \a — 27 ^ m * were calculated in single 

layer system ir*2i. 
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The first two terms are C-S term and Coulomb in- 
teraction term for + gauge field which take exactly the 
same forms as in a single layer system^. The third and 
the fourth terms are non-relativistic Maxwell terms for 
+ and — modes respectively. The last two terms couple 
+ mode to — mode. Integrating out + modes leads to 
e 2 _q 3 (a^) 2 + e-(e-uj 2 + x-Q 2 )(i<l a o a t) wmcn are sub- 
leading to the Maxwell term of clq,cl[. In fact, these 
terms break Time reversal and Parity, in principle, a C-S 
term iqa^a^ will be generated under RG sense. How- 
ever, the coefficient of this generated C-S term could be 
so small that it can be neglected except at experimentally 
unattainable low temperatures. 

(5) Neutral Gapless modes: 

For simplicity, we only consider the balanced case 
and will comment on im-balanced case later. Note that 
V x aT — 2nSp where Sp = 5p2{x) — Spi(x) is the rel- 
ative density fluctuation of the two layers. Introduc- 
ing a variable <f> which is conjugate to Sp(x), namely 
[4>{x), Sp{x')] = iS(x — x 1 ), we can write a spin-wave 
Hamiltonian density: 

n = \x^{5p) 2 + \ P s{Vcp) 2 (9) 

If Sp is treated as a continuous variable, then cj) is a 
free field varying from — oo to oo. By integrating out Sp, 
we get the <f> representation: 

^ = \xs(d T cp) 2 + ^ S (V0) 2 (10) 

Integrating out (f>, we get an effective action density 
in the Sp representation which is dual to the above <j) 
representation: 

C P = \ X ^\5pf + \p-\- q )\Spf (11) 

Plugging the constraint V x a~ = 2nSp into EqnfTTl 
we get: 

£ tt = i(p;V+xrV)K-) 2 as) 

This is consistent with the well-known fact that a 
pure 2 + 1 dimensional U(l) gauge field is dual to a 
2 + 1 dimensional Gaussian model which does not have 
any topological excitations^. Comparing Eqn[12] with 
EqnfS] ( for simplicity, we take Ni = N2 ), we get 
p s = fmj*/Tr,Xs = [2tt 2 (x + So the spin stiff- 

ness scales as the cyclotron gap and the finite charge gap 
of MCF implies finite spin stiffness. In order to compare 
with experimental data in^£, we have to put back h, c, e, e 
and find the spin-wave velocity: 

v 2 = Ps X s= y -^+ — 7 -== 13 

where n is the total density, I is the magnetic length 
and a ~ 1/137 is the fine structure constant. Note that 



the correct expansion of interlayer Coulomb interaction 
is crucial to get the second term. 

By measuring the transport properties at finite tem- 
perature, the authors in 5 found the activation gap Ea ~ 
OAK. By setting Ea = huj* c and plugging the experi- 
mental parameters n — 5.2 x 10 10 cto~ 2 , d/l = 1.61, e = 
12.6 into Eqn[13l we find that the first term is 1.65 x 
10 10 (cm/s) 2 , the second term is 2.54x 10 12 (cm/s) 2 which 
is two orders of magnitude bigger than the first term. Fi- 
nally, we find v ~ 1.59 x 10 e cm/s which is dominated by 
the second term. This value is in good agreement with 
v ~ 1.4 x 10 6 cm /s found ini. 

(6) Topological excitations: 

As discussed in the previous paragraph, at mean field 
theory, there are four kinds of gapped excitations with to- 
tal fractional charges ±e\ — ±vi and ie^ = ±^2 which, 
for example, can be excited by finite temperature close to 
hu*. But their relative charge distributions between the 
two layers are undetermined. In fact, all the possible ex- 
citations can be characterized by their (a + , a - ) charges 
9-)- F° r example, the four kinds of excitations are 
denoted by (±^1, q-), (±1^2 5 3-)- For v\ = 1/2, they re- 
duce to two sets: (±^,g_). If Sp in Eqn[9]is treated as 
a discrete variable, then < 4> < 2ir is an angle variable. 
q_ must be integers 0, ±1, • ■ • . Exchanging a + leads to 
1/r interaction between the four sets of beasts. While 
exchanging a - leads to logarithmic interactions which 
lead to a bound state between two beasts with opposite 
g_ + 1 0. The energy of this bound state with length L is 
E b = A+ + A_ +q%^ + <? 2 hcj* In L/l where A+ and A_ 
are the core energies of QH and QP respectively. 

An important question to ask is what is the gluing con- 
ditions ( or selection rules ) of (q+ , g_ ) for the realizable 
excitations ? Namely, what is spin ( g_ ) and charge 
( q + ) connection ? Two specific questions are: (1) Is 
there a charge neutral vortex excitation with (0,±1) ? 
Being charge neutral, this kind of excitation is a bosonic 
excitation. (2) Is there a charge ±1/2 and spin neu- 
tral excitation (±1/2,0) ? If they do exist, then the 
QP and QH pair with q_ = have the lowest energy 
tkd* = A + + A_. They decouple from a - gauge field, 
therefore interact with each other only with 1 jr interac- 
tion and are asymptotically free even at T — just like 
those in single layer system. Their charges are evenly 
distributed in the two layers ( namely carry fractional 
charges ±1/4 in each separate layer ). They are decon- 
fined ±1/2 excitations which are completely different ex- 
citations from merons which are confined logarithmically. 

Unfortunately, the spin-charge connection is far from 
obvious in this MCF approach. So little can be said about 
these two possibilities. In the composite boson approach 
to be presented in the next section, both interesting pos- 
sibilities are ruled out. 

(7) Extension to other rilling factors: 

Let's briefly discuss ut = 1/2 bilayer system. It is 
well known that this state is described by Halperin's 331 
stated. In this state, the singular gauge transformation 
U — e i S d2x S d2x ' u c<f3Pc,{3)arg(x-x')pp(x) w ^ ere the matrix 
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U = g ^ i 2 j attached two intralayer flux quanta and 

one interlayer flux quantum to electrons to form Entan- 
gled Composite Fermions ( ECF ). At layer 1 ( layer 2 
), the filling factor of ECF is u{ = ( v\ = |j* ). 
Only when the two layers are identical N\ = N 2 , we get 
v l =v 2 = ^ QH states on both layers, therefore vr = 1/2 
system lacks interlayer coherence. 

(8) Comparison with another version of CF approach: 
It is instructive to compare our MCF picture devel- 
oped in this section with the earlier pictures proposed 
inil. When N% = N 2 , the authors in 1 — attached cf> = 2 
flux quanta of layer 2 to electrons in layer 1 or vice versa 
to form interlayer composite fermions so that at mean 
field theory, CF in each layer form a compressible Fermi 
liquid. They conjectured that a~ gauge field fluctua- 
tion mediates an attractive interaction between CF in 
different layers which leads to a ( likely p-wave) pairing 
instability. This pairing between CF in different layers 
opens an energy gap. But no systematic theory is devel- 
oped along this picture. In the MCF picture studied in 
this section, there is a charge gap which is equal to the 
cyclotron gap even at mean field theory which is robust 
against any gauge field fluctuation. 

(9) Interlayer tunneling: 

The interlayer tunneling term is: 

H t = tc\(x)c 2 (x) + h.c. (14) 
Substituting Eqn|4]into above equation leads to: 

JJ t — ttp\{x)e i $ d2x 'l ar 9{ S ~ S ')p2{x')-arg(x'~x)p 1 (x')]^^^ 

+ h.c. (15) 

This is a very awkward equation to deal with. 

The authors ir*2 pointed out that the tunneling process 
of an electron from one layer to the other corresponds to 
an instanton in the 2+1 dimensional compact QED. 
They applied the results of Polyakov on instanton- anti- 
instanton plasma on 2 + 1 dimensional compact QED and 
found the effective action: 

C = ^( d vX) 2 + ^ cos * ( 16 ) 

where g is the U(l) gauge coupling constant given in 
Eqn llOl In the original work of Polyakov, x is a non- 
compact field — oo < x < oo. The compactness of \ was 
forced in an ad hoc way ir*2. Note that the compactness 
of QED has nothing to do with the compactness of X- 

(10) Summary of limited success of MCF theory: 

We use a Mutual Composite Fermion (MCF) picture 
to explain the interlayer coherent incompressible phase 
at d < d c . This MCF is a generalization of Composite 
Fermion (CF) in single layer QH systems to bilayer QH 
system o 33 i 34 ' 35 i 36 . In the mean field picture, MCF in each 
layer fill exactly v* = 1 MCF Landau level. There are 
four kinds of gapped quasi-particles (QP) and quasi-holes 
(QH) with total fractional charges ±e^ = ±i/i,±e 2 = 



±1/2 ( For v\ = i>2 = 1/2, they reduce to two sets ). 
When considering the fluctuations above the v* = 1 MCF 
QH states, we identify the NGM with linear dispersion 
relation u> ~ vk and determine v in terms of experimental 
measurable quantities. When interpreting the activation 
gap Ea found ir£ in terms of the cyclotron gap of the 
v* = 1 MCF Landau level, we calculate v and find its 
value is in good agreement with the value determined 
in^. Tentatively, we intend to classify all the possible 
excitations in terms of a + and a~ charges (q + ,q_). 

(11) Serious problems with MCF approach: 

Despite the success of the MCF approach mentioned 
above, there are many serious drawbacks of this ap- 
proach. We list some of them in the following. 

(a) The determination of the fractional charges is at 
most intuitive. A convincing determination can only be 
firmly established from the CB approach to be discussed 
in the next section. Furthermore, the MCF still carries 
charge 1, while the QP or QH carry charge ±1/2. An 
extension of Murthy-Shankar formalism 36 in SLQH to 
BLQH may be needed to express physics in terms of these 
±l/2QPandQH. 

(b) It is easy to see that the spin wave dispersion in 
Eqnll3l remains linear u> ~ vk even in the d — ► limit. 
This contradicts with the well established fact that in the 
d — > limit, the linear dispersion relation will be replaced 
by quadratic Ferromagnetic spin-wave dispersion relation 
uj ~ k 2 due to the enlarged SU(2) symmetry at d — > 
0. This is because the flux attachment singular gauge 
transformation Eqn[2] breaks SU(2) symmetry at very 
beginning even in the d — > limit. 

(c) The broken symmetry in the ground state is not ob- 
vious without resorting to the (111) wavefunction. The 
physics of exciton pairing can not be captured. The ori- 
gin of the gapless mode is not clear. 

(d) The compactness of the angle 4> in EanfTUlwas put 
in by hand in an ad hoc way. 

(e) The spin-charge connection in (q + , g_ ) can not be 
determined. 

(f) The interlayer tunneling term can not be derived 
in a straight-forward way. See section III-E. 

(g) In the imbalanced case, there are two MCF cy- 
clotron gaps hw* a at mean field theory. However, there 
is only one charge gap in the system. It is not known 
how to reconcile this discrepancy within MCF approach. 

(h) It is not known how to push the MCF theory fur- 
ther to study the very interesting and novel Pseudo-spin 
density wave state to be discussed in section IV. 

In the following, we will show that the alternative CB 
approach not only can achieve all the results in this sec- 
tion, but also can get rid of all these drawbacks natu- 
rally. Most importantly, it can be used to address the 
novel state: Pseudo-spin density wave state in interme- 
diate distances. 
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III. COMPOSITE BOSON APPROACH: 
EXCITONIC SUPERFLUID STATE 

Composite boson approach originated from Girvin and 
Macdonald's off-diagonal long range order—, formulated 
in terms of Chern-Simon Ginsburg- Landau theory 3 ^. It 
has been successfully applied to Laughlin's series v — 
2^— - It has also been applied to BLQBi°.i±2. Unfor- 
tunately, it may not be applied to study Jain's series at 

p 



2sp±l 



with p ^ 1 and v — 1/2 Fermi Liquid system. 
In this section, we applied CSGL theory to study both 
balanced and im-balanced bi-layer QH system. Instead 
of integrating out the charge degree of freedoms which 
was done in all the previous CB approac h 10 ! 11 , we keep 
charge and spin degree of freedoms on the same footing 
and explicitly stress the spin and charge connection. We 
study how the imbalance affects various physical quanti- 
ties such as spin wave velocity, the meron pair distance 
and energy, the critical in-plane field for commensurate- 
incommensurate transition, etc. 
We can rewrite Hi nt in Eqn[T] as: 

Hint = -SpiVSpi + ^8p 2 V5p 2 + Sp 1 V5p 2 



-Sp + V+5p + + -Sp-VSp. 



(17) 



where V = V n = V 22 = ^f,V = V 12 = ^e"« d , V± = 

Vz ^ v and Sp± = Spi ± Sp 2 . 
Performing a singular gauge transformatio n 31 ! 32 : 

(j) a (x) = e *fd 2 x'arg(S-S')p(S') Ca ^ ^ 

where p(x) — c\{x)c\{x) + c\{x)c 2 {x) is the total density 
of the bi-layer system. Note that this transformation 
treats both ci and c 2 on the same footing ( See Fig. 3). 
This is reasonable only when the distance between the 
two layers is sufficiently small. 

It can be shown that 4> a (x) satisfies all the boson com- 
mutation relations. It can be shown that this SGT is the 
same as that in Eqn[2]if we choose tj> = 1/2 and replace 
both pi and p 2 by p = p\ + p 2 . Because the SGT with 
6 = 1/2 in Eqnj2] is not single- valued, the statistics is 
changed from fermion to boson. 

We can transform the Hamiltonian EqnQ]into the La- 
grangian in Coulomb gauge: 



C = 4>l(d T -iao)<t> a +<t>\{% 
1 



(-ifiV + -A{x) - M{x)) 2 



2m 



d x 1 5 p(x)V+{x — x)5p{x') 



1 



+ - / d 2 x'8p_(x)V-{x - x*)8p-{x') - — a (V x a) 



2tt 









c 










FIG. 3: The flux attachment given by Eqn ll8l There is one 
flux quantum penetrating both layers wherever there is one 
electron. It does not matter this electron is in layer 1 or 2. 
One can compare this Fig. with Fig. la. For simplicity, we 
only show three electrons in the top layer and two electrons 
in the bottom layer. 



spin in the absence of Zeeman term, so the Lagrangian 
has a SU(2) pseudo-spin symmetry. The V- term breaks 
the SU(2) symmetry into U(l) symmetry. In the BLQH 
V- > 0, so the system is in the Easy-plane limit. 

By taking full advantage of the easy-plane anisotropy 
shown in Eqnll9[ we can write the two bosons in terms 
of magnitude and phase 



K = \/Pa + Spa,e V 



(20) 



The boson commutation relations imply that 
[5p a (x),6b(x)] = ihS a bS(x - x' ). 

We can see EanlT51 has a local U(l) gauge symmetry 
a -> a + x.a^ -> a M + d^x and also a global U(l) 
symmetry di — ► 0% + x, 9 2 — » 9 2 — \- We denote the 
.^symmetry by U(1) L x U(1) G . 

Absorbing the external gauge potential A into a and 
substituting EqnJ2D]into Eqn[Tni we get: 



(19f = \ d rSp + + \p+id T 6 + + ^6p+(d T d+ - 2a ) 



In Coulomb gauge, integrating out ao leads to the con- 
straint: V x a — 2n(j>\(j> a . Note that if setting V- = 0, 
then the above equation is identical to a single layer with 



2 ' 2 

Pa + Sp a 



2m 



(W a - a) 



i 

2 l 

l -Sp+V + {q)8p+ ^-a (V 



+ \{S P - + p-)d T 6- + l -8p-V-{q)5p- 
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where 5p = 5p\ ± 5p2, p± = Pi ± P2, Q± = 9i±9 2 - They 
satisfy the commutation relations [Sp a (x), 9p(x')\ — 
2ih5 a p5(x — x') where a, (3 — ±. 

It is easy to see that the symmetry is U(l)+ x J7(l)_ 
where U(l)+ is a local gauge symmetry, while U(l)- is 
a global symmetry. We also note that it is Sp~ + p- = 
(dpi — 6p2) + {pi — p-i) = p\ — P2 which is conjugate to the 
phase namely, [Sp~ (x)+p-, 9-(x')] — 2ihS(x—x'). In 
the following, for the simplicity of notation, we redefine 
5p~ = p\—p2, then the last two terms in Eqn l21l becomes: 



only consider the low energy spin-wave excitation. The 
charge sector ( 9 + mode ) and spin sector ( 9~ mode ) 
are essentially decoupled. 

(1 ) Off-diagonal algebraic order in the charge sector: 
The charge sector is essentially the same as the CSGL 
action in BLQH. Using the constraint a t = 27rS q p , ne- 
glecting vortex excitations in the ground state and inte- 
grating out Sp + leads to the effective action of 9 + : 



±6p-V-(q)6 P 



c, 



1 



e+{-q,-u)[ 



•<4 



V-P- =V-(p! 



h z 5p- (22) 
pi) plays a role like a Zee- 



V+(q) 



47T 2 p 1 



-(?» (26) 



where h z ■ 
man field. 

By expressing the spatial gradient term in Eqn l21l in 
terms of (<5p+,#_|_) and (5p- 7 9-), we find 



-%q 2 V + (q) and u c = ^ is the cyclotron 
frequency. 

From Eqn l26l we can find the equal time correlator of 



C = i8p+(~d T e + -a ) + £-[±V6 + + ±(v 1 

+ \5p + V+{q)5p+ - ^a (V x a) 
i lit 



v 2 )V0- - a 



2 < 6 + (-q)6+(q) > 



2^ 



= 2x^ + 0(1) 
q 2 q 



ts p - are - + £L { l ve _ f 



\strV-($)Sf 



where / = <±v\Vi which is equal to 1 at the balanced case. 

Plus the extra terms due to the magnitude fluctuations 
in the spatial gradient terms: 



/i z <5/(23)vhich leads to the algebraic order: 



2m 



Sp^ 



1 



V6L 



(24) 



Note that in Eqn[24]both 8p + and 8p~ couple to terms 
with two spatial gradients, therefore can be dropped rel- 
ative to the terms in Eqn l23l which contains just one tem- 
poral gradient. 

In the following, we will discuss balanced and imbal- 
anced cases separately: 



< 9 + (-q,-Lj)9 + (q,Lj) > 
(27) 

F (28) 



Note that if we define 9+ 



< e 



i(9+(x)-9+(y)) 



> = 



= 0+/2, then 
| a _^ 1 y 2 which takes exactly the 
same form as that in v = 1 SLQH. However, when consid- 
ering vortex excitations to be discussed in the following 
e ^+(x) mav no t De single valued, therefore the fundamen- 
tal angle variable is + instead of 9+ . 
(2) Spin-wave excitation: 

While the spin sector has a neutral gapless mode. In- 
tegrating out 5p~ leads to 



2V-(q) y 2 ' 2m V 2 



(29) 



Balanced case v\ = U2 = 1/2 



where the dispersion relation of spin wave can be ex- 
tracted: 



Putting v\ = V2 = 1/2 and h z = into EqnJ5 
the Lagrangian in the balanced case: 

C = i5p + ( l -d T 9+ - ao ) + ^(iw+-a) 2 
+ ^Sp+V + (q)5p+ - ^a (V x a) 

l Z7T 



-dp d T t 



p 



1. 



-(-V0_ 
2m V 2 



1 



we get 



6p-V-(q)6p-(25) 



In the balanced case, the symmetry is enlarged to 
U(1)l x U(1)g x Z2 where the global Z2 symmetry is 
the exchange symmetry between layer 1 and layer 2. At 
temperatures much lower than the vortex excitation en- 
ergy, we can neglect vortex configurations in Eqn l25l and 



= {^V-(q)]q 2 = v 2 (q)q 2 



In the long wavelength limit: 



V-(q) 



ne 
e 



(d - -d 2 q + -d 3 q 2 + 



2 1 6 
The spin wave velocity is: 



2 2/ ^ P m2 j e2 
v = v f q = 0) = d = — 1 

m e me 



(30) 



qd^.1 (31) 



Eqn l32l shows that the spin wave velocity should in- 
crease as Vd when d < d c . At d = 0, v = 0. This 
is expected, because at d = the U(1)g symmetry is 
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enlarged to SU(2)g, the spin wave of isotropic ferro- 
magnet u> ~ k 2 . By plugging the experimental param- 
eters m ~ 0.07m e which is the band mass of GaAs, 
p = 5.2 x 10 10 cto~ 2 , d/l = 1.61, e = 12.6 into Eqn|52l we 
find that v ~ 1.14 x 10 7 cto/s. This value is about 8 times 
larger than the experimental value. Although Quantum 
fluctuations will renormalize down the spin stiffness from 
Pbare = P to p e // < P, it is known that CSGL theory can 
not give precise numerical values on energy gaps even in 
SLQH. In fact, the spin stiffness p s which is defined as 
4f(V0_) 2 in Eqnf2H] should be determined by the inter- 
layer Coulomb interaction instead of being dependent of 
the band mass to. 

(3) Topological excitations: 

Any topological excitations are characterized by two 
winding numbers A9i = 27rmi, A#2 = 2nm2, or equiva- 



lently, A9 + = 2-n(mi + m-i) 



2ttto+, A6>_ 



27t(toi 



m<i) — 27tto_. It is important to realize that the two fun- 
damental angles are 9\ , 9i instead of 9 +,9-. mi , rcii are 
two independent integers, while to + , to_ are not, because 
to + — m_ = 2to,2 which has to be an even integer. 

There are following 4 kinds of topological excitations: 
A0i = ±2tt, A6»2 = or A6»i = 0, A6> 2 = ±2tt. Namely 
(toi,TO2) = (±1,0) or (mi,TO2) = (0, ±1). They corre- 
spond to inserting one flux quantum in layer 1 or 2, in 
the same or opposite direction as the external magnetic 
field. Let's classify all the topological excitations in terms 
of {q, m_) where charge q is the fractional charge of the 
topological excitations in the following table. 



(toi,to 2 ) 


(1,0) 


(-1,0) 


(0,1) 


(o,-i) 


TO_ 


l 


-1 


-1 


1 


TO_|_ 


l 


-1 


1 


-1 


q 


1/2 


-1/2 


1/2 


-1/2 



Table 1: The fractional charge in the balanced case 

The fractional charges in Table 1 were determined from 
the constraint V x a = 2it5p and the finiteness of the 
energy in the charge sector: 



q 



= — I a ■ dl= — x - I V6»+ • dl= -to+ (33) 
2vr / 2i 2/ 2 + K ' 



to_ = ±1 gives the vorticities which lead to logarith- 
mic interaction between the merons, while q = ±1/2 
lead to 1/r interaction. Therefore merons are confined 
into the following two possible pairs at low temperature. 

(1) Charge neutral pairs: (±1/2, ±1) or (±1/2,^1)- The 
NGM will turn into charge neutral pairs at large wavevec- 
tors ( or short distance ). The pair behaves as a boson. 

(2) Charge 1 pair (1/2, ±1) or charge —1 pair (—1/2, ±1). 
The pair behaves as a fermion. They are the lowest 
charged excitations in BLQH and the main dissipation 
sources for the charge transports. A duality transforma- 
tion ( See EqnJUl) can be easily performed to express low 
energy physics in terms of the dynamics of these topolog- 
ical excitations. There is a possible Kosterlitz-Thouless 
(KT) transition above which the meron pairs are liber- 
ated into free meron. 



The MCF picture in the last section points to two in- 
teresting possibilities (1) There maybe Charge neutral 
bosonic excitations with (0,±1): Note that m + = im- 
plies toi = — 7?i 2 = to and to_ = 2to. For to = 1, it 
corresponds to inserting one flux quantum in layer 1 in 
one direction and one flux quantum in layer 2 in the op- 
posite direction. So only (0, ±2) exist, while (0, ±1) do 
not exist. (2) There may be deconfined ( or free ) 1/2 
charged excitations. Because any excitations with non- 
vanishing to_ will be confined, so any deconfined excita- 
tions must have m_ = which implies mi — to-2 = to 
and to + = 2to. Wc find the charge q = \m + = to must 
be an integer. This proof rigorously rules out the pos- 
sibility of the existence of deconfined fractional charges. 
We conclude that any deconfined charge must be an in- 
tegral charge, mi = TO2 = 1 corresponds to inserting one 
flux quantum through both layers which is conventional 
charge 1 excitation. Splitting the two fluxes will turn 
into a meron pair with the same charge. 

(4) Comparison with the LLL+HF approach 
It is constructive to compare the spin ( — ) sector of 
Eqnf2S] with the EPQHF Hamiltonian achieved by the 
microscopic LLL±HF approach ini 2 -: 



■ P 

'2 
PA 



£ = i'-A(m) ■ d T m + j3 m {m z ) 2 - Cqm z (-q)m z (q) 



where Vrs 



termined by 

J °° dxx 2 exp( 



A 
the 



m, [3 m ~ 
2 exp(-x 2 /2) 



C 



(34) 



167re ' 

3 de- 



16V27reZ 

intralayer interaction and pe = 
-x 2 /2 — dx/l) is determined by the 
Note that all these numbers are 



16-n-tl 

interlayer interaction 
achieved by assuming that the ground state wavefunc- 
tion is Halpcrin's 111 wavefunction even at finite d. How- 
ever as shown i n 39 ' 65 , the wavefunction at any finite d is 
qualitatively different from Halperin's (111) wavefunction 
which is good only at d — 0. So the numbers calculated 
by the LLL±HF based on (111) wavefunctions may not 
even have qualitatively correct distance dependence. 

In the above equation, the first term is the Berry phase 
term, the second term is the mass ( or capacitance ) term, 
this term leads to the easy-plane anisotropy which sup- 
presses the relative density fluctuations between the two 
layers, the third term is nonanalytic in the wave vector 
due to the long range nature of the Coulomb interaction, 
the fourth term is the spin stiffness term for m z and the 
fifth term is the spin stiffness term for easy-plane. At 

2 

d = , P m = C = Q,p A = Pe = 16 3^ el ; tnen £ in 
Eqnj34| reduces to the SU(2) symmetric QH ferromag- 
net as it should be. Note that the value of pA = pe at 
d = is exact, because the ground state wavefunction 
is exactly the Halperin (111) wavefunction, while at any 
finite d, the estimates of pa ^ Pe by HF approximation 
may be crude, because we still do not know what is the 
exact groundstate wavefunction which may be qualita- 
tively different from the (111) wavefunctio n 12 ! 39 . In the 
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presence of the easy-plane ( (3 m ) term, C and pa terms 
are subleading, therefore, can be dropped in the long 
wavelength limit. However, they are still very important 
if there is a instability happening at finite wavevector as 
shown in section III-B. 

If taking the symmetry breaking direction to be along 
the x direction, we can write ra x = yl — m 2 cos 8, m y = 
iy/1 — m 2 sin#, m z with m z ~ 0,8 ~ 0. Substituting the 
parameterizations into the Berry phase term in Eqnl34| 
we find the Berry phase term to be ijm z d T 8, if identi- 
fying m z = 26 p-, it is the same as the linear derivative 
( ^6p-d T 8~ ) term in the spin sector of Eqn[25l How- 
ever, we can see the leading term in Eqnl31l which leads 
to the capacitive term is ~ d, but in Eqn(33J it is ~ d 2 , 
while the Monte-Carlo simulation in spherical geometry 
ir*2£ indicates it is ~ d. The subleading term in Eqnl3~H 
is non-analytic ~ q instead of ~ q 2 , this is due to the 
long-range behavior of the Coulomb interaction. This 
non-analytic term is the same as that in Eqn l34l if iden- 
tifying m z = 25 p- . The third term in Eqn l31l leads to 
a (Vm z ) 2 term with a coefficient ~ d 3 , while the coeffi- 
cient of (Vm 2 ) 2 term in Eqnl34l is pa which approaches 

2 

the constant 16 ^n^ el as d — ► as dictated by SU(2) sym- 
metry at d = 0. The difficulty to recover SU(2) sym- 
metric limit at d = from the CB approach is due to 
that the decomposition EqnJ20]in our CB approach takes 
advantage of the easy-plane anisotropy from the very be- 
ginning. This is similar to Abelian bosonization versus 
Non-Abelian bosonization in one dimensional Luttinger 
liquid or multi-channel Kondo model^. The last ( Pe ) 
term in EcmlMI becomes ^(V6») 2 . 

By this detailed comparison between the CB approach 
and the microscopic LLL+HF approach, we find that the 
two approaches lead to exactly the same functional form 
in the spin sector, some coefficients such as the Berry 
phase term and the non- analytic term are the same, 
while some other coefficients such as the easy-plane term 
and spin stiffness term are not. Unfortunately it is very 
difficult to incorporate the LLL projection in the CB ap- 
proach. As suggested mi 2 -, we should simply take some 
coefficients in CB approach as phenomenological values 
to be fitted into microscopic LLL+HF calculations or 
numerical calculations or eventually experimental data. 
The advantage of CB approach is that it also keeps the 
charge + sector explicitly, therefore treat the QH effects 
in the charge sector and the inter-layer phase coherence 
in the spin sector at the same footing. While the charge 
sector in the LLL+HF approach is completely integrated 
out. 

Most of the results achieved in this subsection were 
achieved before ^ 11 ' 12 in microscopic calculations where 
the charge fluctuations were integrated out, the LLL pro- 
jections were explicitly performed and HF approxima- 
tions were made. Here, we reproduce these old results in 
a very simple way which keeps both spin and charge in 
the same footing and bring out the spin-charge connec- 
tion in a very transparent way. We also classify all the 



possible excitations in this effective CB approach. 

In the next section, we will look at the effects of im- 
balance. 



B. Im-balanced case v\ 7^ vi 

(1 ) Off-diagonal algebraic order and Spin-wave excita- 
tion: 

In the im-balanced case, the second term in Eqn[23] in- 
cludes the coupling between spin sector and charge sector 
even when neglecting vortex excitations. Expanding this 
term, we find the effective action of the coupled 8 + and 
8- modes: 



-8 + (-q,-Uj)[ 



V+(q) 

LQ 2 
V-{q) 



iir 2 p 1 



-(<?>) 



-q 2 ]8-(q,uj) 



P ( 
I^ 1 



V2)q 2 8-{-q, -w)6+(q,w) 



(35) 



where we safely dropped a linear derivative term in (9_ in 
Eqn. [23] in the Interlayer Coherent QH state. However, 
the linear derivative term will be shown to play important 
role in the PSDW to be discussed in section IV. B. From 
Eqn l351 we can identify the three propagators: 



< 



< 



< 



> 



> = 



(|^)(c 2 + ^yy 2 

lu a + lu 2 {uj 2 + v 2 q 2 ) + fv 2 q 2 uj 2 

W-(q)(L, 2 +L0 2 ) 



uj 4 + uJ 2 {uJq + v 2 q 2 ) + fv 2 q 2 UJ 2 

-Myi - v 2 )V~{q)u 2 q 
lu a + w 2 (u> 2 + v 2 q 2 ) + fv 2 q 2 uj 2 



(36) 



where / = 4^i^ 2 < 1, w| = uj 2 + ^q V+(q), the cyclotron 
frequency lo c = and the spin wave velocity in the bal- 
anced case v 2 — v 2 {q) — ^V-(q) were defined in EqnlHOl 
Performing the frequency integral of the first equation 
in Eqn l36l carefully, we find the leading term of the equal 
time correlator of 8+ stays the same as the balanced case 
EqnETl 

2tt 1 

< 8 + {~q)8 + {q) >= 2 x — + O(-) (37) 
T 1 

which leads to the same algebraic order exponent 2 as in 
the balanced case Eqnl2"51 

We conclude that the algebraic order in the charge 
sector is independent of the imbalance. This maybe ex- 
pected, because the total filling factor z/y = 1 stays the 
same. 

In the q,Lu — > limit, we can extract the leading terms 
of the 8-8- propagator: 



< 



>= 



4^-(g) 
2 + fv 2 q 2 



(38) 
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where v 2 — v 2 (q = 0) and we can identify the spin wave 
velocity in the im-balanced case: 



fv% = 4viV 2 Vq =4^i(1 



(39) 



which shows that the spin- wave velocity attains its maxi- 
mum at the balanced case and decreases parabolically as 
the im-balance increases. The corresponding KT transi- 
tion temperature Txt also decreases parabolically as the 
im-balance increases. 

Similarly, in the q, to — > limit, we can extract the 
leading terms in the 8+8- propagator: 



> 



<m- >= ~ 4F - ( f ;r/ 2) = - (*-*)<* 

u> z + jv q 

(40) 

which shows that the behavior of < 8+8- > is dictated 
by that of < 8-8- > instead of < 8 + 8 + >. 

When the vortex excitations to be discussed in the fol- 
lowing are included, the spin wave velocity will be renor- 
malized down. As analyzed in detail in the last subsec- 
tion, it is hard to incorporate the Lowest Landau Level 
(LLL) projection in the CB approach, so the spin wave 
velocity can only taken as a phenomenological parameter 
to be fitted to the microscopic LLL calculations or exper- 
iments, its precise dependence on imbalance can only be 
determined by experiments. 

(2) Topological excitations 

There are following 4 kinds of topological excitations: 
A0i = ±2tt, A8 2 = or A0i = 0, A8 2 = ±2tt. Namely 
(mi, ma) = (±1,0) or (mi, ma) = (0,±1). We can clas- 
sify all the possible topological excitations in terms of 
(g, m_) in the following table. 



(mi, ma) 


(1,0) 


(-1,0) 


(0,1) 


(0,-1) 


m,- 


1 


-1 


-1 


1 


m+ 


1 


-1 


1 


-1 


q 


V\ 


-v\ 


V2 


-vi 



Tabic 2: The fractional charge in im-balanced case 

The fractional charges in table 2 were determined from 
the constraint V x a = 2n8p and the finiteness of the 
energy in the charge sector: 

q = -!- la-dl= -!- x \ /[V0+ + [y x - i*)V0~] ■ dl 



2 TV 

1 

2 



2tt 2 

m+ + (v\ - i/ 2 )m. 



(41) 



In contrast to the balanced case where q only de- 
pends on m+ , q depends on m+ , rri- and the filling fac- 
tors ^1,^2- Just like in balance case, any deconfincd 
excitations with (m_ = 0, m+ = 2m) have charges 
q = hm+ = m which must be integers. While the charge 
of (m+ = 0, TTi— = 2m) is q = {yi — V2)m which is charge 
neutral only at the balanced case. 

The merons listed in table 2 are confined into the fol- 
lowing two possible pairs at low temperature. (1) Charge 
neutral pairs: (±^i,±l) or (±^a, ±1)- They behave as 



bosons. The NGM will turn into charge neutral pairs at 
large wavevectors. (2) Charge 1 pair (fi,l) + [yi,— 1) 
or charge —1 pair (— v\,— 1) + (— fa,l). They behave 
as fermions. They are the lowest charged excitations in 
BLQH and the main dissipation sources for the charge 
transports. 

In appendix A, we will calculate the meron fractional 
charges from its trial wavefunction and find they are the 
same as those listed in table 2. 

(3) Dual action 

A duality transformation can be easily performed to 
express low energy physics in terms of the dynamics 
of these topological excitations. Performing the du- 
ality transformation on EqnJ55] leads to the dual ac- 
tion in terms of the vortex degree of freedoms = 
■^tiiv\d il d v 6± = J^ 1 ± J^ 2 and the corresponding dual 
gauge fields &*: 

C d = -inb+e^xdvb^ - iAf u e^ x d„b+ + iirb+J v u 



777 I — * — * 

+ - dnbtf + -(V x b+)V + (q)(V X b+) 



i^ s ^u\d u b x + iirb^ - h z (\7 x b ) 
^(9 tt 6 - - d Q b-) 2 + i(V x b-)V-(q)(W x b~) 
jjiyx - V2)(d b- ~ d b-){d (3 b+ - 8 b+) (42) 



where Af^ = A\ ± A 2 ^ are the two source fields. The 
+ sector stands for the charge sector which is essentially 
the same as SLQH. While the spin ( or — ) sector of this 
dual action takes similar form as a 3D superconductor in 
an external magnetic field h z shown in Eqn l63l The last 
term which stands for the coupling between the charge 
and the spin sector can be shown to be irrelevant in the 
ILCQH state. The spin wave velocity in Eqnj39]can also 
be easily extracted from this dual action. This dual ac- 
tion was used to derive the ground state, quasi-hole and 
a pair of quasi-hole wavefunctions ir*^. 

It is constructive to compare this dual action derived 
from the CB approach with the action derived from MCF 
approach Eqn[8l We find the following three differences: 
(1) The topological vortex degree of freedoms are 
missing in Eqn[5J These vortex degree of freedoms are 
needed to make the variable 4> in Eqn[TU] to be an angle 
variable < < 2tt. (2) The x+>X- terms in Eqn[8] 
are extra spurious terms. These extra spurious terms 
break SU{2) symmetry even in the d — * limit. (3) 
The linear term —h z (\7 x b~ ) is missing in EqnJS] This 
linear term is not important in the interlayer coherent 
excitonic superfluid state, but it is very important in the 
in-coherent excitonic solid state to be discussed in the 
following subsection. Indeed, if we drop all the x+iX- 
terms, use bare mass and also add the topological vortex 
currents and the linear term by hands in Eqn|8l then 
EqnJS] will be identical to Eqn. 221 We conclude that 
Eqn. [23] and [15] from CB approach are the correct and 
complete effective actions. 
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(4 ) Energy of a meron pair 

We can also look at the static energy of a charge 1 
meron pair consisting a meron with charge v\ and charge 
v 2 separated by a distance R: 

E mp = E lc + E 2c +'^^+2nf Ps0 ln^- (43) 
it He 

where Ei c ,E 2c are the core energies of meron 1 ( charge 
v\ ) and meron 2 ( charge v 2 ) respectively, R c is the core 
size of an isolated meron, p S Q = ^ is the spin stiffness 
at the balanced case. 

Minimizing E mp with respect to R leads to an optimal 

2 

separation: R Q — - which is independent of the im- 
balance. Namely, the optimal separation of a meron pair 
remains the same as one tunes the imbalance. Plugging 
R into Eqn[43] leads to the optimal energy of a meron 
pair: 

e 2 R 
E = E lc + E 2c + z/i)(— + 87rp s0 In -£) (44) 

H H c 

Because of logarithmic dependence on R c , we can ne- 
glect the v\ dependence in R c , then the second term in 
Eqn. f53] decreases parabolically as im-balance increases. 
Unfortunately, it is hard to know the v\ dependence of 
E\ c ,E 2c and R c from an effective theory, so the depen- 
dence of the energy gap E on V\ is still unknown. 



D. Summary 

In summary of this section, the effective CB theory 
can be used to lead to correct algebraic off-diagonal long 
range order, low energy functional forms, classification of 
both spin wave and topological excitations. But it may 
not be used to get correct values of many physical mea- 
surable quantities such as spin wave velocity, its precise 
dependence on the imbalance and the energy of a meron 
pair. The EPQFM approach is a microscopic one which 
takes care of LLL projection from the very beginning. 
However, the charge sector was explicitly projected out, 
the connection and coupling between the charge sector 
which displays Fractional Quantum Hall effect and the 
spin sector which displays interlayer phase coherence was 
not obvious in this approach. While in the CB theory, it 
is hard to incorporate the LLL projection ( see however 
the attempt made ir>2& ), some parameters can only be 
taken as phenomenological parameters to be fitted into 
the microscopic LLL+HF calculations or experimental 
data as discussed in section III. The two approaches are 
complimentary to each other. The biggest advantage of 
the CB approach is that it can be extended to capture 
competing orders even at microscopic length scales, so 
can also be used to describe novel phases and phase tran- 
sitions as the layer distance changes to be discussed in 
the next section. 



C. Interlayer tunneling 

Note that in contrast to MCF, the same phase factor 
appears in the singular gauge transformation Eqn ll8l for 
the two layers a — 1,2 which get canceled exactly in H t : 

H t = tcj)\{x)4> 2 {x) + h.c. = 2ty/pip 2 cos(#i - 2 ) 
= 2tp^viu 2 cos 0_ (45) 

where the exciton operator is c{c2 = 0i^2 = 
pTjviV 2 e ie - ~ 

In the presence of in-plane magnetic field Bn = 
(B x , By), the effective Lagrangian is: 

+ 2tpy/viv 2 cos{9- - Q a x a ) (46) 

where x a = (x,y),Q a = (-^, 

In balanced case, it was found that when the ap- 
plied in-plane magnetic field is larger than a critical field 
B > Sn ~ (to/pso) 1 / 2 , there is a phase transition from 
a commensurate state to an incommensurate state with 
broken translational symmetry. When B > BT,, there 
is a finite temperature KT transition which restores the 
translation symmetry by means of dislocations in the do- 
main wall structure in the incommensurate phase. 

As can be seen from EqndHl Ps ~ v\v 2 , while t ~ 
y^iT^, so the critical field Bn — (tj p s ) x l 2 ~ {viv^ 1 ^ 
increases as one tunes the im-balance. 



IV. COMPOSITE BOSON APPROACH: 
PSEUDO-SPIN DENSITY WAVE STATE 

The discussions in the last section are in the homoge- 
neous exciton superfluid state, now we will study the in- 
stability of this ESF state as the distance is increased. As 
stated in the introduction, when the distance is smaller 
than d c \, the system is in the ESF state, while it is suf- 
ficiently large, the system becomes two weakly coupled 
v = 1/2 Composite Fermion Fermi Liquid (FL) state. 
There could be a direct first order transition between 
the two states. However, the experimental observations 
that both zero voltage tunneling peak 4 ' and the Hall drag 
resistivity^ develop very gradually when d ~ d c \ suggest 
the transition at d — d c \ is a 2nd order phase transi- 
tion. Although there are very little dissipations in both 
the ESF and FL, the experiment^ discovered strong en- 
hancement of drag and dissipations in an intermediate 
distance regime. These experimental observations sug- 
gest that there must be intermediate phases between the 
two phases. In this section, we will study the nature of 
the intermediate phase. 

In this section, by combining the composite boson 
method developed in the last section with the Ginzburg- 
Landau theory developed in2&, we explore the similarity 
between the 4 He system and the BLQH step by step. In 
the dual density representation of Eqn l21l we propose a 
Pseudo-spin density wave ( PSDW ) state which breaks 
translational symmetry as the candidate of the intermedi- 
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ate phase. We also construct a novel effective Ginzburg- 
Landau theory to study the novel transition from the 
ESF to the PSDW. Inspired by recent possible discovery 
of supersolid in 4 i? o 23 ' 24 i 28 , we also discuss important 
physical consequences of quantum fluctuation generated 
vacancies in the PSDW state. Just like in the last sec- 
tion, we study the balanced case first, then the effects of 
imbalance. 



Balanced case 



(1 ) Collapsing of magneto-roton minimum 
As shown in III-A-4, the spin stiffness and the 
V-(q) in Eqn[2H] should be replaced by the effective 
ones calculated by the LLL+HF approximation: p E and 
V E {q). 



2V E (q) [ 2° T 



1. 



(47) 



where Ve(q) — a — bq + cq 2 , a,b,c > 0, a ~ d 2 ,6 ~ d 2 , but 



j 2 , a, b, c > 0, a ~ d 2 , b ~ d 2 , 
c remains a constant at small distances^ 2 -. As explained 
in the last subsection, the non-analytic term is due to the 
long-range Coulomb interaction which survives the LLL 
projection. 

In order to study the instability of the ESF to the 
formation of PSDW, one need to write the dispersion 
relation Eqn l30l to include higher orders of momentum: 



= [2p E V E (q)]q 2 = q 2 (a-bq + cq 2 ) 



(48) 



As shown ir*2&, the advantage to extend the dispersion 
relation beyond the leading order is that the QGL action 
can even capture possible phase transitions between com- 
peting orders due to competing interactions on micro- 
scopic length scales. As explained in the last section, in 
the LLL-HF approach iiJ* 2 ., these coefficients were found 
to be a ~ d 2 , b ~ d 2 , but c remains a constant at small 



distances. By looking at the two conditions ui 
and 



,o=0 

^-\ q=qo = 0, it can be shown that the disper- 
sion curve Eqn05] indeed has the shape shown in Fig. 5a. 
When b ~ d 2 < b c = 2 v / ac ~ d, the magneto-roton min- 
imum has a gap at q = qo = \Ja/c ~ d, the system in 
the ESF state, this is always the case when the distance 
d is sufficiently small ( Fig. 4). However, when b = b c , 
the magneto-roton minimum collapses at q = qo which 
signifies the instability of the ESF to some density wave 
formation as shown in2&. When b ~ d 2 > b c = ~ d, 

the minimum drops to negative and is replaced by the 
stable density wave formation, the system gets to the 
PSDW state, this is always the case when the distance d 
is sufficiently large ( Fig. 4). 

The phenomenon of the collapsing of the magento- 
roton minimum as the distance increases was also seen in 
the numerical calculations in the LLL+HF approaches 56 . 
It was estimated that qol ~ 1, so the lattice constant of 




b~cT 



2 (ac) ~ d 



dci 



FIG. 4: The critical distance d a \. 
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FIG. 5: The zero temperature phase diagram in the balanced 
case as the distance between the two layers increases. ESF 
where < ip 0, < uq >= stands for excitonic superfluid, 
PSDW where < ip >= 0, < >j= stands for pseudo- 
spin density wave phase, FL stands for Fermi Liquid, (a) 
Energy dispersion relation U)(q) in these phases, (b) Vs(g) 
in these phases. The cross in the PSDW means the negative 
minimum value of Ve (?) is replaced by the PSDW. The order 
parameters are also shown. In fact, the instability happens 
before the minimum touches zero. 



the PSDW is of the same order of magnetic length I which 
is ~ lOOA. The critical distance d c i is also of the same 
order of the magnetic length. 

(2) Instability in the density channel 

From Eqns l25l29l30l we can see that it is the original 



instability in V E {q) 



bq + cq which leads to the 



magneto-roton minimum in the Fig. 5. By looking at the 
two conditions V E (q)\ q=qo = and \ q = qo = 0, it 

is easy to see that V E (q) indeed has the shape shown in 
Fig. 5b. When b ~ d 2 < b c = 2y/ac ~ d, the minimum of 
V E (q) at q = qo = \J ajc ~ d has a gap, the system is in 
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the ESF state, this is always the case when the distance d 
is sufficiently small. However, when b = b c , the minimum 
collapses and S(q) diverges at q = q ( Fig. 5a), which 
signifies the instability of the ESF to an exciton normal 
solid (PSDW) formation. When b ~ d 2 > b c = 2yfa~c ~ 
d, the minimum drops to negative, the system gets to the 
PSDW state, this is always the case when the distance 
d is sufficiently large. These conditions are essentially 
the same as those achieved by looking at the dispersion 
relation ui in Fig. 5a. 

(3) Effective action in the dual density representation 
In Eqnf2H Sp- was integrated out in favor of the phase 
field 9- . Because the original instability comes from the 
density-density interaction, it is important to do the op- 
posite: integrating out the phase field in favor of the 
density operator in the original action EqnfJHl Neglect- 
ing the vortex excitations in and integrating out the 
6- in Eqng7| leads to: 

C[Sp-] = bp-(-q,-UJ n )[-^ + V E (q)]5p-(q,u n ) 
2 2p E q l 

(49) 

where we can identify the dynamic pseudo-spin density- 
density correlation function: 



weak interactions among the pseudo-spin density 5p- : 



S-(q,oj n ) =< Sp (-q,-w n )6p (q,u n ) >'- 



2p E q 2 



where v 2 (q) = 2p E V E (q) is the spin wave velocity defined 
in EqnEDl 

From the pole of the dynamic density-density correla- 
tion function, we can identify the speed of sound wave 
which is exactly the same as the spin wave velocity. This 
should not be too surprising. As shown in liquid 4 He, 
the speed of sound is exactly the same as the phonon 
velocity. Here, in the context of excitonic superfluid, we 
explicitly prove that the sound speed is indeed the same 
as the spin wave velocity. 

From the analytical continuation iuo n — > cu + id in 
Eqn l50[ we can identify the dynamic structure factor: 
S-(q,uj) = S-(q)5(uj - v(q)q) where S-(q) = p E qir/v(q) 
is the equal time pseudo-spin correlation function shown 
in Fig. 5a. As q — > 0, S-(q) — ► q. The Feymann rela- 
tion in BLQH which relates the dispersion relation to 
the equal time structure factor is 



u(q) 



PE^q 

S-(q) 



(51) 



which takes exactly the same form as the Feymann 
relation in superfluid 4 He. Obviously, the V E (q) in 
the Fig. 5b leads to the magneto-roton dispersion uj 2 = 
q 2 V E (q) in the Fig.5a. 

Because the instability happens near q = qo instead of 
q = 0, so the transition in Fig5 is not driven by vortex un- 
binding transitions like in 3D XY model, so the vortices 
remain tightly bound near the transition^ 5 .. So integrat- 
ing out the vortex excitations in will only generate 



2~ r v " " ,l 2p E q 2 
+ u{5p_f + w(6p-) 6 + ■■■ 



V E {q))Sp (q,u) 

(52) 



where the momentum and frequency conservation in the 
quartic and sixth order term is assumed. 

In sharp contrast to the conventional classical normal 
liquid (NL) to normal solid (NS) transition^, the possi- 
ble cubic interaction term (5p-) 3 is forbidden by the Zi 
exchange symmetry between the two layers Sp- — > —Sp- . 
Note that the (w n /q) 2 term in the first term stands for 
the quantum fluctuations of 5p- which is absent in the 
classical NL to NS transition. The density representa- 
tion Eqnf52]is dual to the phase representation Eqnl29l 
However, the phase representation Eqn l29l contains ex- 
plicitly the superfluid order parameter ip_ ~ e l ®~ which 
can be used to characterize the superfluid order in the 
ESF phase. All the results in the ESF state are achieved 
from phase representation in the last section. While in 
Eqn l521 the only signature of the superfluid phonon mode 
is encoded in the density sound mode, because the order 
parameter is integrated out, the superfluid order is 
hidden, so it is not as powerful as the phase represen- 
tation in describing the ESF state. However, as shown 
in the following, when describing the transition from the 
ESF to the PSDW, the density representation Eqn[53lhas 
a big advantage over the phase representation. 

Expanding V E (q) near the minimum qo in the Fig. 3 
leads Eqn l52l to the quantunm Ginsburg-Landau action 
to describe the ESF to the PSDW transition: 



C[5 P -] = l5p-[A p u>l + r + ctf-qtf]5p- 
+ u(5p-) 4 + w{6p_f + --- 



(53) 



where A n ~ ^ ? which is non-critical across the tran- 

sition. 

Because near q , the coefficient A p of u 2 in Eqnl53lis 
a constant in the density representation, while the cor- 
responding quantity Ag ~ (q) in the phase represen- 
tation Eqn|27| is divergent, so Eqn[57] breaks down as 
q — > q$ and may not be used to describe the ESF to the 
PSDW transition. 

It is constructive to compare the QGL Eqn[53] to de- 
scribe the ESF to the PSDW in 2 + 1 dimensional BLQH 
with that in 3 + 1 dimensional 4 He to describe the su- 
perfluid to normal solid transition^: 



C[Sn] = -6n[A 

- w(8n) 3 



'l+r + c{q 2 
u(8n) 4 + --- 



ql) 2 ]8n 



(54) 



which is non-critical 

Psqg 



where r ~ p cl — p and A p 
across the transition. 

The most crucial difference between Eqnj53] and 
EqnJSHis the absence of the cubic term in Eqnf53]and the 
presence of the cubic term in Eqn l54l As shown in the 
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FIG. 6: 2,4,6 shortest reciprocal lattice vectors in a (a) stripe 
(b) square lattice (c) triangular lattice 

following, this crucial difference make the lowest energy 
lattice structure different in the two systems. Drawing 
the analogy of the familiar He4 phase diagram at T = 
discussed in2£, we may assume the divergence of S~ (q) at 
q = qo leads to the PSDW state in the Fig. 5. In Eqnl53l 
r is the gap of Ve(<?) at the minimum which tunes the 
transition from the ESF to the PSDW. This is a Bra- 
zovskii type transition^ described by a n = 1 component 
(d+ l,d) quantum Lifshitz action ( with d = 2 in BLQH 
Due to the absence of the cubic term, the transition 
is 2nd order at mean field level. However, as shown by 
renormalization group analysis in2S, the fluctuations will 
drive the transition into a 1st order transition. 
(4 ) Lattice structures of the PSDW phase 
In the ESF, r > and < <5p_ >= is uniform. In the 

PSDW, r < and < 5p- >= Eg «(G)e 4c5 s , ti(0) = 
takes a lattice structure. In two dimension, the most 
common lattices are Id stripe embedded in two dimen- 
sion, square and triangular lattices. The shortest recip- 
rocal lattice vectors of the three lattices were shown in 
Fig. 5 (a) and (b) and (c) respectively. 

(a) P=2: Qi = —Q2 = Q are a pair of anti-nodal 
points ( Fig. 6a ). They are the two shortest reciprocal 
lattice vectors generating a 1 dimensional lattice embed- 
ded in a 2 dimensional system. The Id stripe will lead to 
transport anisotropy which was not seen in experiments, 
so will not be considered any more in the following. 

(b) P=4: Q 3 = -Qi, Qi = -Q 2 , Qi ■ Q-z = 0, Q t ,i = 
1, 2, 3, 4 form the 4 corners of a square ( Fig. 6b ). They 
are the four shortest reciprocal lattice vectors generating 
a 2 dimensional square lattice. 

(c) P=6. Qi,i — 1,2,3,4,5,6 form the 6 corners of a 
hexagon ( Fig. 6c ). They consist of the 6 shortest recipro- 
cal lattice vectors generating a 2 dimensional triangular 
lattice. 

If thinking Eqnf53] as at 2 + 1 dimension, it can be 
shown that the triangular lattice ( Fig. 6c ) is the favorite 
lattice for the normal solid. As shown in the appendix B, 
due to the cubic term, the triangular lattice is also the 
favorite lattice for the Wegner crystal (WC) in the QH to 
WC transition in the SLQH. However, due to the absence 
of the cubic term in Eqn l531 in the following, we will show 
that the favorite lattice is the square lattice. At mean 
field level, we can ignore the lo dependence in Eqn l53l 
Substituting < 5p_ >= Ti(C?)e iG£ into EqnrSTil leads 



to: 

G G 
G 

where u could be positive or negative, v > is always 
positive to keep the system stable. 

From Eqnl55[ we can compare the ground state energy 
of the two most commonly seen lattices: square lattice ( 
Fig. 6b ) and triangular lattice ( Fig. 6c ). Square ( trian- 
gular ) lattice has m = 4(m = 6) shortest reciprocal 
lattice vectors G — 2n/a where a is the lattice constant 
in a given layer ( Fig. 6). For the quadratic and quartic 
terms, all the contributions come from the paired recip- 
rocal lattice vectors ( Fig. 6b and 6c ). For the sixth 
order term, in square lattice, all the contributions are 
still from the paired reciprocal lattice vectors ( Fig. 6b 
), however, in triangular lattice, there is an additional 
contribution from a triangle where G\ + G2 + G3 = ( 
Fig. 6c ). It is this extra contribution which make a cru- 
cial difference between the square lattice and the trian- 
gular lattice. Following^, after scaling tiq — > mr 1 / 2 ng 
to make the quadratic term the same for both lattices, 
then EqnfSniis simplified to: 

fa = \r 6 \n d \ 2 + u a \n d \ A + v a \n e \ 6 (56) 

where u\j = ma = 3u for both square and triangular 
lattice and v\j = 18u,wa = 18u + 20/3v for square and 
triangular lattice respectively. Obviously . The mean 
field phase diagram of EqnJ55] is well known: If it > 
( u < ), there is 2nd ( 1st ) order transition, there is 
a tri-critical point at u = 0. Minimizing / with respect 
to riQ leads to ng — n a = (— 2u a + y 1 ^u 2 a — 6v a t)/6v a 
which holds for both u a > and u a < 0. Obvi- 
ously, 7ja 7^ Tin. From Eqn l56l we can see that be- 
cause v/\ > Vrj, for any given n, /□(«■) < /a(ti), then 
fa( n n) < /p(tia) < /a («a), namely, the square lattice 
is the favorite lattice . So it has three elastic constants 
instead of two. Neglecting zero-point quantum fluctua- 
tions, < 5p- >= J2i 8(x — Ri) — — Ri — where 
the I = 5(01 + 02) is the shift of the square lattice in the 
bottom layer relative to that in the top layer ( Fig. 2). 

This PSDW state not only breaks the translational 
symmetry, but also the Z2 exchange symmetry. It is 
very rare to get a 2d square lattice, because it is not a 
close packed lattice. Due to the special Z2 symmetry of 
BLQH, we show it indeed can be realized in BLQH. The 
system is compressible with gapless phonon excitations 
determined by the 3 elastic constants. 

(5) Light scattering from the PSDW state 
It is very interesting to see if soft X-ray or light scat- 
tering experiments^ can directly test the existence of the 
PSDW when d c \ < d < d c2 - The in-plane light scatter- 
ing intensity from the PSDW is I{K) ~ \S{K)\ 2 where 
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(a) (b) 

FIG. 7: (a) The charge distribution of the PSDW in a square 
lattice. The "up " pseudo-spins take sublattice A, while the 
"down" pseudo-spins take sublattice B. (b) In general, there 
are quantum and thermal fluctuation generated vacancies in 
both layers denoted by a O- 

S(K) is the structure factor and K = —(i+J) is the 2d 
in-plane reciprocal lattice vector of the square PSDW. 
For f = §(?+/), S(K) = 1 - e iJ? ' r = 1 - (-l) m +« which 
is and 2 for even and odd m + n respectively. The in- 
plane light scattering experiments 00 probe the intralayer 
shift I in the Fig. 7, in principle, it can be performed in 
the BLQH. Note that the large zero-point quantum fluc- 
tuations may diminish the light scattering intensity by a 
Debye- Waller factor—. 

The out of-plane light scattering intensity from the 
PSDW follows the Bragg diffraction law 2<isin6' = nX 
where 9 is the angle between the light ray and the 2d 
plane and the A is the wavelength of the light. It doubles 
or vanishes when n is an integer or a half integer. The 
out of-plane light scattering experiments probe the inter- 
layer distance d in the Fig. 7, but it is quite difficult to 
perform due to several GaAs/AlGaAs layers above the 
2d electron gas. 

(6) Vacancies, disorders and Coulomb Drag in the 
PSDW state 

In principle, the 8p+ mode in Eqn l23l should also be in- 
cluded. It stands for the translational (or sliding ) motion 
of the PSDW lattice. However, any weak disorders will 
pin this PSDW state to make it an insulating state (INS). 
The insulating state also has a charge gap which has a 
completely different origin than that in the QH state. Be- 
cause there are charge gaps on both sides of the ESF/QH 
to PSDW/INS transition, it is valid to integrate out the 
charge 5p+ sector and focus on the spin sector in studying 
the two phases and the transition in Fig. 5. Disorders will 
smear the 1st order transition from the ESF to the PSDW 
into a 2nd order transition. It was argued in-i that dis- 
orders in real samples are so strong that they may even 
have destroyed the ESF state, so they may also trans- 
fer the long range lattice orders of the PSDW into short 
range ones. This fact makes the observation of the lattice 
structure by light scattering difficult. Being an insulat- 
ing state, the PSDW state will not show any quantized 
Hall plateau and any zero-bias interlayer tunneling peak. 
The two square lattices in the top and bottom layer are 



locked together, so it will show huge Coulomb drag. How- 
ever, vacancies generated by the large zero-point quan- 
tum fluctuations may play important roles in the drag 
( Fig. 2b ). As the distance increases to the critical dis- 
tance d c \ in Fig.l, the ESF turns into the PSDW whose 
lattice constant a = V^LtZ is completely fixed by the fill- 
ing factor v\ = 1/2, so it may not completely match the 
instability point l/q®. Due to this slight mismatch, the 
resulting PSDW is likely to be an in- commensurate solid 
where the total number of sites N s may not be the same 
as the total number of sites N even at T = 0. As the 
distance increases further d c \ < d < d C 2 in Fig.l, the 
PSDW lattice constant is still locked at a = ^/AttI. As- 
suming zero-point quantum fluctuations favor vacancies 
over interstitials, so there are vacancies no even at T = 
in each layer. At finite temperature, there are also ther- 
mally generated vacancies n a (T) ~ e~ Al '/ T where A„ is 
the vacancy excitation energy. So the total number of 
vacancies at any T is n v (T) = no + n a (T). Obviously, 
the vacancies in top layer are strongly correlated with 
those in the bottom layer. The drag in the weakly cou- 
pled FL regime d > d C 2 was shown to be very small and 
goes like T 4 / 3 at low Ti 3 .. The drag in the ESF regime 
d < d c i is also small and goes like e~ A<3ff / T at low T. 
The drags in both states are due to momentum trans- 
fer between electrons in the two layers. However, in the 
PSDW, the holes are hopping on the square lattice, as 
shown in 4 ^, the microscopic origin of the drag is due to 
the correlated character of hopping transports of the va- 
cancies between the active and passive layers, it leads to 
a very large drag. We can estimate the resistance in the 
active layer as R ~ l/n v (T). As shown inii, the drag 
resistance in the passive layer is Rd ~ aoR where ocd 
should not be too small. This temperature dependence 
is indeed consistent with that found in the experiment 20 : 
starting from 2§0mK, Rd increases exponentially until 
to bOmK and then saturates. This behavior is marked 
different than that at both small and large distance where 
the system is in the ESF and FL regimes respectively. We 
suggest that in the presence of disorders, all the proper- 
ties of the PSDW are consistent with the experimental 
observations ini^^ on the intermediate phase. The ef- 
fects of very small imbalance on the PSDW phase will be 
investigated in IV-B-3 and was found to naturally explain 
the small imbalance experiment in££. All these facts sug- 
gest that the intermediate phase can be identified as the 
PSDW. 

(7) Finite temperature phase diagram 

Due to the spontaneous translational symmetry break- 
ing in the PSDW, there are also Goldstone modes in the 
spin sector which are the two lattice phonons u(r) in 
2 dimension. The translational order is characterized by 
n Q ffi = e l ® m ' u , the orientational order is characterized 
by ipor — e l4t ^ for the square lattice where <f> = \z ■ V x u 
( in contrast, ip or — e 1 ^ for the hexagonal lattice ). It is 
known that at any finite temperature at d — 2, there is 
no true long range order except possible long-range orien- 
tational order. The square lattice has 3 elastic constants 
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FIG. 8: The Finite temperature phase diagram in the 
balanced case as the distance between the two layers in- 
creases. ESF stands for excitonic superfluid with only al- 
gebraic ODLRO. PSDW stands for Pseudo-spin density wave 
with only algebraic translational order and true orientational 
order. The EH stands for the exciton hexatic phase with 
only algebraic orientational order. FL stands for Fermi Liq- 
uid. Tkt is the KT transition temperature driven by vortex 
unbinding in the phase 9-. T v is the vector Coulomb gas 
transition temperature driven by dislocation unbinding. T 3 
is the scalar Coulomb gas transition temperature driven by 
disclination unbinding. 



( due to higher symmetry, the hexagonal lattice only 
has 2 ). We assume the Kosterlitz-Thouless-Halperin- 
Nelson- Young ( KTHNY ) theoryii designed for hexago- 
nal lattice also work for the PSDW square lattice in Fig. 7. 
When interlayer distance d = 0, the symmetry is SU(2), 
the broken symmetry state at T = in Fig. 3 is imme- 
diately destroyed at any finite T in Fig. 8. The system 
becomes the exciton liquid (EL). When < d < d c i, the 
true long-range ordered PSDW in Fig. 5 has only alge- 
braic ODLRO in Fig. 8, there is a finite temperature KT 
transition above which the system is in the EL. When 
d c i < d < d C 2, the true long-range ordered translational 
order in Fig. 5 has only algebraic one, but true orienta- 
tional order. There could also be a dislocation driven 
transition at T = T v in the universality class of 2d vec- 
tor Coulomb gas with the correlation length exponent 
v v = 0.37. The system is in the Exciton Hexatic ( EH 
)phase with short range translational order and algebraic 
orientational order. There could also be a disclination 
driven transition at T = T a in the universality class of 
KT transition ( scalar Coulomb gas ) with the correlation 
length exponent v s = 1/2. The system becomes the EL 
phase. The finite temperature phase diagram is shown 
in Fig.8. 

(8) The PSDW to the FL transition 

In Fig. 5, as the distance increases to d C 2, the PSDW 
will melt into the two weakly coupled FL. The flux at- 
tachment transformation EqnJTS] which treats the elec- 
trons in the two layers on the same footing breaks down 
and a completely different flux attachment transforma- 
tion within a single layer to transform an electron into a 



Composite Fermion Eqn[5]must be used^. At d < d C 2, 
it is the composite boson to feel zero field, they either 
condense into the phase ordered ESF or the density or- 
dered PSDW states, while at d > d C 2, it is the compos- 
ite fermion to feel zero field, they form a Fermi surface. 
Unfortunately, because completely different actions are 
needed in the two sides of the transition at d = d C 2 , this 
quantum phase transition can not be addressed in this 
paper and remains an interesting open question. 

(9) General arguments on the existence of the PSDW 
state 

When the distance between the two layers is very large, 
intralayer correlations dominate, a crystalline state is not 
favorable. When the distance is sufficiently small, due 
to the strong interlayer correlation, electrons in the top 
layer pair with the holes in the bottom layer to form 
excitons, it is essentially a bosonic system in the pseudo- 
spin channel. As shown in all the experiments^^, the 
strong interlayer correlations go well beyond the regime 
of the ESF phase. At T = 0, any bosonic system has 
to be ordered in some way, namely, breaks some sort of 
symmetry. If it is a homogeneous state, then it is phases 
ordered which breaks some internal symmetry such as 
global U(l) symmetry, this is the ESF phase where the 
kinetic energy of the excitons dominate. It happens at 
very small distance d < d c ±. If it is an inhomogeneous 
state, then it is density ordered which break translational 
symmetry, this is the PSDW phase where the potential 
energy of the excitons dominate. This happens at some 
intermediate distances d c \ < d < d C 2- There is no other 
possibilities. Some exotic state such as exciton liquid 
which does not break any symmetry is very unlikely^ 

In the following, we will discuss the effects of the im- 
balance on Fig. 5 



B. Imbalanced case 

(1 ) The shape of the magnetoroton minimum 

In order to see how the imbalance affects the shape 
of the dispersion curve shown in the Fig. 5a, we need to 
push Eqn[35]to the order q 4 : 

u;f mb = f^-f(l-f)^- (57) 

where / = ^v\V2 and u> 2 is defined in Eqn l481 

Scaling out the common factor /, one can see a and 
b stay the same, while c is reduced by the imbalance. 
Therefore qo increases, namely, the lattice constant a of 
the PSDW may decrease ( See Fig. 10 ). While b c de- 
creases, namely, the critical distance d c ± also decreases. 
Unfortunately, as warned in Sec. III. D, the CB approach 
may not describe the magnetoroton minimum shape de- 
pendence on imbalance correctly. In the subsection (3), 
we will approach the phase boundary from the PSDW 
side. 

(2) Dual density representation in the ESF state 
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Just like in the balanced case, in order to focus on the 
original instability, we need to keep 8p-. The form of 
Eqn l35l before integrating out Sp- is 



1 



LO 2 + LO 2 . 



£c = ^9 + (-q,-Lj){ : - \_ - ]0+(q,w) 

+ jz.i 1 '! ~ v2)q 2 0-(-q, -u)d+{q,u>) 

+ \5p-d T e- + ^(^V0_) 2 + ±5p-V-tf)5p- - K 

Obviously, integrating out 5p- and dropping out the 
linear derivative term in 9~ , we recover Eqn l35l 

By defining 5p^ — ^- + 5p'_, we can absorb the last lin- 
ear term in Eqnl58l into the quadratic term. Integrating 
out both + and leads to: 



2 2 i 2 

2 pq 2 io 2 + fw* 

where we can identify the dynamic pseudo-spin density- 
density correlation function S'(q,u>) =< 5p'_5p'_ > : 



< Sp'_Sp'_ >- 



uj 4 + lu 2 (lu 2 + v 2 q 2 ) + fv 2 q 2 u> 2 



(60) 



where v 2 = v 2 (q) was defined in Eqn[30l Obviously 
S'(q,uj = 0) = v ^. 

Again, From the pole of the dynamic density-density 
correlation function Eqn l601 we can identify the speed of 
sound wave in the imbalanced case which is exactly the 
same as the spin wave velocity in the imbalanced case 
Eqn l57l achieved from the phase representation. From 
both the phase representation in Sec. Ill and the density 
representation in this subsection, we conclude that im- 
balance is irrelevant in the ESF side, but it may play 
important role in the PSDW side to be discussed in the 
following. 

(3) The quantum phase transitions driven by small im- 
balances in the PSDW state 

As shown in section III, the imbalance is irrelevant in 
the ESF side, but we expect it is important in the PSDW 
side. Starting from the ESF side, it would be useful to 
calculate how the mageto-roton minimum depends on the 
imbalance, this was attempted in (1). As said in (1), 
the precise dependence of d c \ on the imbalance is hard 
to achieve from the CB effective theory. Here we take 
a different strategy: starting from the PSDW side and 
studying how a small imbalance affects the PSDW. If the 
imbalance is sufficiently small, we expect the C-PSDW 
at the balanced case in Fig. 7a is a very good reference 
state. In the following, we ignore the small number of 
vacancies in the PSDW shown in Fig. 7b. Because it is 
a square lattice with the "up" pseudospins taking sub- 
lattice A and the "down" pseudospins taking sublattice 
B, it is reasonable to start from a lattice model from the 
PSDW side. If we think the PSDW as a charge density 
wave (CDW ) of bosons at half filling on a square lattice, 



then we can view the ESF to the PSDW as a superfluid 
to CDW transition in a boson Hubbard model of hard 
core bosons near half filling hopping on square lattice of 
bosons: 



H = -t ^2 (b\bj + h.c.) - fi^rii 

<ij> i 

+ Vx fa - 1/2) ( nj - 1/2) 

<ij> 

+ V 2 fa-l/2)(n fe -l/2) + --- (61) 



«ik» 

where S+ = &t = c \c 2 ,S- = b = c\c u S z = \(c{c x - 

c\c2) = 6^6 — 1/2 are the pseudo-spin density and bo- 
son operators. At the total filling factor v<r = 1, we can 
impose the local constraint c\c\ + c 2 C2 = 1. Hj = b\bi 
is the boson density, t is the nearest neighbor hopping 
amplitude, V±, V2 are the nearest and next nearest neigh- 
bor repulsive interactions between the bosons. The • ■ ■ 
may include further neighbor interactions. Because of the 
long-rang Coulomb interaction in Eqn l231 it is important 
to keep all the long-range interactions in the lattice model 
EqiifrSJ If the chemical potential p — 0, the bosons are 
at the half filling < m >= 1/2 which corresponds to the 
balanced case v — 1/2. The particle- hole symmetry of 
Eqnll9l corresponds to the Z 2 exchange symmetry of the 
BLQH. If the chemical potential p ^ 0, the bosons are 
slightly away from the half filling which corresponds to 
the slightly imbalanced case. 

The boson Hubbard model Eqn. [HI] in square lattice 
at generic commensurate filling factors / = p/q ( p,q 
are relative prime numbers ) were systematically studied 
ir*££ by performing the charge-vortex duality transfor- 
mation. Recently, we applied the dual approach to study 
extended boson Hubbard model Eqn[6T] in bipartite lat- 
tices such as square and honeycomb lattice which maybe 
directly relevant to the Helium 4 and Hydrogen adsorbed 
on various structures of substrate a 52 i 53 . The key result 
achieved in this paper is that there are two consecu- 
tive transitions at zero temperature driven by the cov- 
erage : a Commensurate-Charge Density Wave (CDW) 
at half filling to a narrow window of supersolid, then to an 
Incommensurate-CDW. In the Ising limit, the supersolid 
is a CDW supersolid; whereas in the easy-plane limit, 
it is a valence bond supersolid. The first transition is a 
second order one in the same universality as the Mott 
insulator to the superfluid transition, therefore has the 
exact critical exponents z = 2, v = 1/2, r) = ( with 
possible logarithmic corrections ). The second is a first 
order transition. In the following, we simply straightfor- 
wardly apply the results near half filling in the square 
lattice achieved in^i to the present problem on square 
lattice near half filling. At q = 2, there are two dual 
vortex fields ip a and ?/>e,. Moving slightly away from half 
filling / = 1/2 corresponds to adding a small mean dual 
magnetic field H ~ Sf = f — 1/2 in the dual action. The 
most general action invariant under all the MSG trans- 
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formation laws upto quartic terms i a 50 ' 51 : 



-^ty.v\d v A\ ~ 2nSfS fJ , T 



2^2 



Oi—a/b 

+ 7o(|^a| 2 + |^| 2 ) 2 "7i(l^| 2 -|^| Z ) 

where A^ is a non-compact £7(1) gauge field. Upto the 
quartic level, Eqn l62l is the same in square lattice and 
in honeycomb lattice. If r > 0, the system is in the 
superfiuid state < i\>i >= for every I = a/b. If r < 0, 
the system is in the insulating state < ipi >^ for at 
least one I. We assume r < in Eqnl6"2l so the system 
is in the insulating state. 71 > ( 71 < ) corresponds 
to the Ising ( or Easy-plane ) limit. The insulating state 
takes the CDW state ( or valence bond solid (VBS) state 
)■ 

In the balanced case Sf = 0, the SF to the VBS 
transition in the easy plane limit was argued to be 
2nd order through a novel deconfined quantum critical 
point 54 . However, the boson Hubbard model EanlTOl on 
the PSDW side corresponds to the Ising limit in the dual 
model EqnEl therefore 71 > 0. The SF to the CDW 
transition in the Ising limit is first order. This is consis- 
tent with the first order ESF to PSDW transition driven 
by the collapsing of magnetoroton minimum studied in 
Sec.IV-A. 

In the CDW order side, the mean field solution is ip a — 
1, ipb = or vice versa. In the slightly imbalance case 
Sf 7^ 0, setting ipb = in Eqn[52] leads to: 



C = \(d^-iA^ a \ 2 +r\^ a \ 2 + lQ \ij a \ 4 
+ ~(V x A-2vr5/) 2 



(63) 



where pa — Wa4>a should be interpreted as the vacancy 
number, while the vortices in its phase winding are inter- 
preted as boson number. Of course, a negative imbalance 
can simply achieved by a particle hole transformation 
^ ^t, */->-*/ ill EqnJBSJ 

Eqn l63l has the structure identical to the conventional 
q = 1 component Ginzburg-Landau model for a type II 
" superconductor " in a "magnetic" field. It was well 
known that as the magnetic field increases, there arc two 
first order phase transitions: H < H c \, the system is in 
the Messiner phase, H c \ < H < H C 2, it is in the vortex 
lattice phase, H > H C 2 it is in the normal phase. In the 
present boson problem with the nearest neighbor inter- 
action V\ > in EqnlBTl which stabilizes the CDW state 
at / = 1/2, this corresponds to C-CDW to IC-CDW to 
superfiuid transition shown in Fig. 9a. Transferring back 
to the original BLQH problem, the small imbalance will 
first drive the C-PSDW to the In-commensurate normal 
solid (IC-PSDW), then drive a 1st order transition from 
the IC-PSDW to the ESF shown in Fig.9b. 

As shown in Fig. 9b, the bias voltage increases, the im- 
balance will first introduce interstitials in the top layer 
and vacancies in the bottom layer, namely, turn the C- 
PSDW into the IC-PSDW whose charge distributions are 
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FIG. 9: (a) The phase diagram of the boson Hubbard model 
Eqn ll9l slightly away from the half filling, /i is the chemical 
potential, (b) The bias voltage V versus distance d phase di- 
agram at zero temperature. IC-PSDW stands for the incom- 
mensurate PSDW. The dashed line is the experimental path 
investigated All the transitions are first order transi- 

tions. 




(a) 



(b) 



FIG. 10: (a) The charge distribution of the PSDW in a square 
lattice. The dashed line stand for "up" pseudo-spins which 
take sublattice A, while "down" pseudo-spins take sublattice 
B. (b) The charge distribution of the IC-PSDW in a square 
lattice. The dotted line stand for "inverted" pseudo-spins 
which play the role of interstitials in the top layer or vacan- 
cies in the bottom layer. The number of interstitials in the 
top layer is equal to that of vacancies in the bottom layer. 
Compare with Fig. 7 



shown in Fig. 7 and 10 respectively, then the whole IC- 
PSDW melts into the ESF through a 1st order transition. 

The dashed line in Fig. 9b was investigated in a re- 
cent experiment 58 . But the first phase transition in 
Fig. 9b was not paid attention in the experiment where 
the phase diagram was drawn against fixed charge im- 
balance h z instead of fixed bias voltage V. So the C- 
PSDW phase was crushed into the horizontal axis. A 
simple mean field argument leads to the linear scaling 
of the second transition line V ~ d — d c \. A parabolic 
behavior h? z ~ d — d c \ was found for the shape of the sec- 
ond transition at very small imbalances. Disorders will 
smear all the 1st order transitions in Fig. 9b into 2nd or- 
der transitions. The disorders may also transfer the long 
range lattice orders of the C-PSDW and IC-PSDW into 
short range ones. The fact make the observation of the 
commensurate and incommensurate lattice structures by 
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light scattering experiment 60 difficult. We expect the dis- 
orders may transform the linear behavior V ~ d — d c i in 
clean case to the parabolic one h 2 z ~ d— d c \ in the dirty 
case. In the presence of disorders, all the properties of 
the C-PSDW and IC-PSDW are consistent with all the 
experimental observations in*£ on the intermediate phase 
at small imbalance. 

When the distance of the two layers is further increased 
to larger than a second critical distance d C 2, then all the 
signature of the inter layer coherent state are lost, the 
two layers are decoupled into two separate v = 1/2 CF 
Fermi liquid state ( Fig. lb ). We expect that there is 
a level crossing and associated first order transition at 
d C 2- When d > d C 2, increasing the bias voltage may not 
transform the two decoupled FL state back into the ESF. 

(4 ) Large imbalance 

At large imbalance, all the states in Fig. 5 at the bal- 
anced case are not good reference states anymore. Fig. 
7 is valid only for small imbalance such that v\ , v<2 do 
not fall onto any fractional Quantum Hall plateaus for 
separate layers. For large imbalance such as two layers 
with V\ = 1/3,^2 = 2/3 still with vt — 1, then the two 
weakly coupled layers show v\ = 1/3 and — 2/3 frac- 
tional Quantum Hall states separately when d > d C 2- 
Because there are intralayer gaps, they are more ro- 
bust against the interlayer correlations, we expect d C 2 
for v\ = 1/3, V2 = 2/3 to be smaller than v\ — V2 = 1/2. 
If we bring them closer, there could be a direct 1st order 
transition at d c \ = d C 2 from the weakly coupled pair to 
the ESF, or through a narrow regime of IC-PSDW. The 
extreme unbalanced case V\ -C 1, zaj = 1 — v% ~ 1 where 
coupled electron and hole Wigner crystal (EHWC) was 
discussed in^. As shown in the appendix B, we believe 
this EHWC should take a triangular lattice in contrast 
to the square lattice at the balanced case. 



C. Comparison with earlier work 

Although the ESF phase and FL phase at the two ex- 
treme distances are well established, the picture of how 
the ESF phase evolves into the two weakly-coupled FL 
states is still not clear, namely, the nature of the interme- 
diate phase at d c \ < d < d C 2 is still under debate. There 
are previous work also based on the instability due to 
magnetoroton minimum collapsing 69 . These work pro- 
posed different kinds of translational symmetry break- 
ing ground states as candidates of the intermediate state 
. All these previous work used either HF approxima- 
tion or trial wavefunctions approximation. Especially, 
Cote, Brey and Macdonald in 6 - 9 - found that the lowest 
energy lattice structure of the PSDW is square lattice. 
The microscopic calculations are powerful and accurate 
in studing quantum Hall states with a bulk gap in SLQH 
or BLQH systems, but its accuracy and power drop con- 
siderably in the v = 1/2 FL system in SLQH and the 
v = 1 BLQH system with a broken symmetry state and 
its associated Goldstone mode. The Quantum Ginsburg- 



Landau theory presented in this paper is complementary 
to and goes well beyond the previous microscopic calcula- 
tions. As shown in the text, our QGL theory Eqn l53l verv 
quickly and firmly leads to the conclusion that the square 
lattice is the most favorite lattice. We also explicitly 
point out the hopping mechanism for the large Coulomb 
drag in the PSDW state. The two approaches are com- 
plementary to each other and reach similar conclusions. 
Furthermore, the QGL theory circumvents the difficulty 
associated with the unknown wavefunction at any finite 
dM and treat both the interlayer and the intralayer cor- 
relations on the same footing, so can be used to capture 
competing orders on microscopic length scales and natu- 
rally leads to the PSDW as the intermediate state which 
breaks translational symmetry. The theory puts the ESF 
state and the PSDW state on the same footing, character- 
ize the symmetry breaking patterns in the two states by 
corresponding order parameters and describe the univer- 
sality class of the quantum phase transition between the 
two states. It can also be used to determine the nature, 
properties and lattice structure of the PSDW state. The 
properties of the PSDW proposed in this paper are also 
consistent with earlier experimental observations in2°-. 

By the microscopic LLL+ HF approach, the authors 
i n 55 i 56 i 57 claimed that claimed that the transition at 
h z = 0, d = d c i is an instability through a 1st order 
transition to a pseudospin density wave state driven by 
the gap closing of magneto-roton minimum at a finite 
wave-vector. Starting from the ESF side, their numeri- 
cal results indicated that the imbalance increases the spin 
stiffness and also the critical distance d c \. All the calcu- 
lations in LLL+HF approach assumes that the ground 
state wavefunction at any finite d is still the (111) wave- 
function. However as shown i n 39 i 65 , the wavefunction 
at any finite d is qualitatively different from the (111) 
wavefunction which is good only at d = 0. So the excita- 
tion spectra in the ESF side calculated by the LLL+HF 
based on the (111) wavefunction may not have quantita- 
tively correct distance dependence. The LLL+HF theory 
inS 6 . completely breaks down in the PSDW side. Further- 
more, in this LLL+HF approach, the charge fluctuations 
are completely integrated out, therefore can not address 
the interplay between the QH effects in the charge sec- 
tor and the interlayer phase coherence in the spin sector. 
The CB field theory approach in this paper circumvents 
this difficulty associated with the unknown wavefunction 
at any finite d (see Sec.V) and was used to show how the 
pseudo-spin density wave state is formed and its proper- 
ties. We look at the effects of small imbalance from the 
PSDW side and map its effect as a chemical potential 
of a hard core bosons hopping on a square lattice near 
half filling, namely, mapping Fig. 9a to Fig. 9b. Our up- 
per phase boundary in Fig. 9b is consistent with that in 56 . 
We also worked out the lattice structures of the C-PSDW 
and IC-PSDW and the whole physical picture along the 
dashed line in Fig. 9b. However, as pointed out in Sec. Ill, 
it is hard to incorporate the LLL projection in the CB. 
Some parameters can only be estimated in the LLL+HF 
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approaches. So the CB approach in this paper and the 
microscopic LLL+HF approach in££ are complementary 
to each other. 

There are other proposals which are not based on 
the collapsing of magnetoroton minimum. Recently, a 
trial wavefunction involving the coexistence of composite 
bosons and composite fermions at finite distance was pro- 
posed in§£. However, it seems to the author that there 
is no phase transitions in this CF-CB coexistence trial 
wavefunction. This is in-consistent with the experimen- 
tal observations. As further elucidated in Sec.V wave- 
function approach is unlikely to be efficient in describ- 
ing transitions. The authors ir* 6 - 7 - argued that the state 
d c i < d < d c2 is a phase separated state between the 
ESF and the FL. Based on this phase separation picture, 
Wang investigated the effect of the interlayer tunneling 68 . 



V. CONCLUSIONS 

(1 ) Comments on different approaches in SLQH 
Three common approaches to SLQH systems are the 
wavefunction ( or first quantization ) 33 ' 59 , Composite 
Fermion Field Theory (CFFT) 3 !^^ anc } Composite Bo- 
son Field Theory (CBFT) approaches 3 -^. The Wave- 
function approach has been very successfully applied to 
study SLQH at Jain's series at v — 2s p±i ■ CFFT has 
been used to study the CF Fermi liquid at v = 1/2 which 
is at the end point of Jain's series in the limit p — > oo. 
However, there are many problems with CFFT approach 
whose equivalence to the wavefunction approach is still 
not obvious even in SLQH. These problems have been 
vigorously addressed in 3 ^. CBFT approach has so far 
been only limited to Laughlin's series v = -^i+l which is 
only a p = 1 subset of Jain's series. It can not be used in 
any simple way to study p > 1 Jain series and v = 1/2 
CF Fermi liquid. Despite having played very important 
roles historically, except describing topological proper- 
ties elegantly, both CFFT and CBFT have not been able 
to give much new information beyond those achieved in 
the wavefunction approach in describing bulk Quantum 
Hall States. One of main reasons for the success of the 
wavefunction approach is that there is a gap in the bulk, 
suitable wavefunction o 33 i 59 can describe both the ground- 
state and low energy excitations quite accurately. Its ac- 
curacy can be checked easily by exact diagonalization in 
a finite size system whose size is beyond only a few mag- 
netic length. Spherical geometry can be used to get rid 
of edge state effects quite efficiently . However, there are 
at least three exceptions even in SLQH which can not be 
described in the framework of wavefunction approach (1) 
The CF Fermi liquid at v — 1/2. Because it is a gapless 
system, the wavefunction may not be very useful, unfor- 
tunately, a simple field theory description is still lacking. 
(2) The edge state: because the gap vanishes on the edge, 
field theory is very convenient to study the gapless edge 
excitations 7 * (3) The quantum phase transition from QH 
to QH or QH to insulating state. It is known that wave- 



function approach is not convenient to describe any phase 
transition. A Field Theory is the most powerful method 
to study phase transitions. Unfortunately, so far, no field 
theory has been able to describe the transitions except 
in some artificial lattice model a 76 i 77 . In appendix B, we 
will construct a quantum Ginsburg Landau theory to de- 
scribe QH to WC transition. 

(2) Summary of our results on BLQH 

As shown i n 39 i 65 , because of the gapless nature, at any 
finite d, the wavefunction is qualitatively different from 
the (111) wavefunction, so the wavefunction approach 
to BLQH is far less powerful in BLQH than in SLQH. 
So Field theory approaches are much more powerful in 
BLQH than in SLQH, especially in the pseudo-spin sec- 
tor which contains the new phenomena not displayed in 
SLQH. In this paper, we used both CFFT and CBFT 
approach to study the balanced and im-balanced BLQH 
systems. In Section II, we used a MCF approach to study 
balanced and im-balanced BLQH systems. We achieved 
some limited success, but also run into many trouble- 
some problems. We explicitly identified these problems 
and motivated the alternative CB approach. Extension 
of Murthy-Shankar's formalism 36 in SLQH to BLQH is 
not expected to fix these problems. In Section III, we 
developed a simple and effective CB theory to study the 
BLQH. The CB theory naturally fixed all the problems 
suffered in the MCF theory presented in Sec. II. By using 
this CB theory in its phase representation, we first stud- 
ied the excitonic superfluid state and re-derived many 
previous results in a simple and transparent way. The 
theory puts spin and charge degree freedoms in the same 
footing, explicitly brought out the spin-charge connec- 
tion in a straightforward way and classified all the pos- 
sible excitations in a systematic way. We made detailed 
comparisons between the spin sector of the CB theory 
with EPQFM derived from microscopic LLL approach. 
Although some parameters in the spin sector can only 
be taken as phenomenological parameters, the functional 
form is identical to the EPQFM. In Sec. Ill, using this 
CB theory in its dual density representation, we then an- 
alyzed the instability in the pseudo-spin sector and found 
the magneto-roton minimum collapsing leads to a new 
ground state which breaks the translational symmetry: 
Pseudo-spin density wave ( PSDW ) state. We proposed 
that for balanced BLQH, there are two critical distances 
dd < d C 2- When < d < d c \, the system is in the ESF 
state, when d c i < d < d C 2, the system is in the PSDW 
state, there is a first order transition at d c \ driven by the 
collapsing of magneto-roton minimum in the pseudo-spin 
channel. When d C 2 < d < oo, the PSDW melts into two 
weakly coupled v = 1/2 Composite Fermion Fermi Liq- 
uid (CFFL) state. There is also a first order transition 
at d = d C 2- However, disorders could smear the two first 
order transitions into two second order transitions. The 
transition from the ESF to the PSDW is unusual because 
the two states break two completely different symmetries: 
the global internal U(X) symmetry and the translational 
symmetry respectively. We construct an effective theory 
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in the dual density representation to describe this novel 
quantum phase transition. The effective action is a n = 1 
component (d+1, d) with d = 2 quantum Lifshitz action. 
The most favorable lattice is a square lattice instead of 
a hexagonal lattice. Because a 2d square lattice is not a 
close packed lattice, so it is hard to be realized experi- 
mentally. This is probably the only experimental realiza- 
tion of a 2d square lattice in any 2 dimensional system. 
The ESF to PSDW transition is similar to the superfhhd 
to normal solid transition in 4 He system with the dis- 
tance playing the role of the pressure. The correlated 
hopping of vacancies in the active and passive layers in 
the PSDW state leads to very large and temperature de- 
pendent drag consistent with the experimental data. The 
PSDW could be the true ground state responsible for all 
the experimental observations in intermediate distances. 
We also applied the CB theory to study the effects of 
imbalance on both the ESF and the PSDW side in corre- 
sponding sect. Ill and IV. On the ESF side, as we tune 
the im-balance, the system supports continuously chang- 
ing fractional charges, We also derived the dual action of 
the CB theory in terms of topological currents and dual 
gauge fields. By comparing this dual action with the ac- 
tion derived from MCF theory, we can explicitly identify 
the missing and the arti-facts of the MCF approach in 
section II. On the PSDW side, we map the square lat- 
tice PSDW at the balanced case into a hard core bosons 
hopping on a square lattice at half filling, then adding a 
small imbalance in the BLQH corresponds to adding a 
small chemical potential in the boson model. Through 
this mapping, we find the imbalance drives two quan- 
tum phase transitions: the first one is from the commen- 
surate PSDW state to an incommensurate PSDW (IC- 
PSDW), the second one is from IC-PSDW state to the 
ESF state. Both transitions are first order transitions. 
We discuss the effects of disorders and finite tempera- 
ture. We compare our results with the previous results 
achieved from the LLL+HF approach, the other propos- 
als on the possible intermediate phases and available ex- 
perimental data. We concluded that CBFT approach 
is superior than CFFT approach in the BLQH having 
ground states with different kinds of broken symmetries. 
It would be interesting to see if the PSDW can also be 
achieved in Tri-layer quantum Hall systems studied inSl 
Only the interlayer coherent (111) state was discussed in 
this paper, it can be easily generalized to other interlayer 
coherent (to, to, to) ( with to odd ) states at total filling 
factors i>t = 1/m. 

(3) General remarks on bosonic and fcrmionic ap- 
proaches 

It is general true that when there is an ordered state 
with broken symmetry and an associated order param- 
eter, bosonic approach is superior than fermionic ap- 
proach. It is well known that the bosonic approach can be 
easily applied to capture competing orders at long wave- 
length scales. As explicitly and systematically shown 
ir*2£ and this paper, it can also be used to capture com- 
peting orders at microscopic length scales. For exam- 



ple, in Quantum Antiferromagnet, fermionic approach 
can only address the disordered phased, while bosonic 
approach^ can address not only the disordered phase, 
but also the ordered phase with broken symmetry and 
the quantum phase transition between the disordered 
phase and ordered phase. Similarly, in quantum spin 
glass, bosonic approach can study spin liquid, spin glass 
and the quantum phase transition between the two, while 
the fermionic representation can only study the spin liq- 
uid phased. In SLQH system, there is only algebraic 
long range order, but no broken symmetry and no true 
ordered state, so CF approach could be very successful. 
Fermionic theory can be even used to describe QH to 
insulating transitions. For example, the QH to insu- 
lating transitions in a periodic potential was described 
in terms of 2 + 1 dimensional Dirac fermion^L 77 . For 
general BLQH (m,m',n) states with to, to' are odd and 
mm' — n 2 ^ , because there are no broken symme- 
try in the ground states and no associated gapless Gold- 
stone modes, we expect the Composite Fermion approach 
works better. For example, (3,3,1) state at vt = 1/2 
can be described in terms of the Entangled Composite 
Fermion (ECF) discussed at section II. While (to, m, to) 
BLQH system has a true U(l) broken symmetry ground 
state and an associated order parameter and a Gold- 
stone mode, the CB approach becomes more effective 
as demonstrated in this paper. Most importantly, we 
use the CB approach to explore the PSDW state which 
breaks the translational symmetry instead of the U(l) 
symmetry. 

(4) The synthesis of roton collapsing in d — 1, 2, 3 

As shown ir>££, due to the long-range Coulomb interac- 
tions between electrons, in the effective low energy theory 
describing the edge reconstruction in the FQHE, there 
are also two low energy sectors at k — and k — k r . It 
is the magneto-roton minimum collapsing at k = k r is 
responsible for edge reconstruction in the edge state of 
FQHE. In one dimensional edge, the roton manifold at 
k = k r becomes two disconnected points. Of course, 1 
dimension is always special. Higher dimensions could 
be completely different. This paper showed that the 
magneto-roton minimum collapsing at d = 2 leads to 
the PSDW state. At 2d, the roton manifold at k — k r is 
a circle. In2&, it was shown that the roton minimum col- 
lapsing driven by pressure could lead to a normal solid 
or a supersolid state in 4 He. At 3d, the roton mani- 
fold at k = k r is a sphere. Combining all these results, 
we find that the roton minimum collapsing could lead 
to novel physics in all possible experimental accessible 
dimensions. 
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APPENDIX A: MERON FRACTIONAL 
CHARGES IN IMBALANCED CASE 

In this appendix, we will evaluate the fractional 
charges of merons in imbalanced case from their wave- 
functions and find they are indeed the same as those 
listed in Table 2. We will first discuss the ground state 
wavefunction, then the meron wavefunction. The bal- 
anced case was discussed in^E. Here we extend the anal- 
ysis to the imbalanced case. 



1. Ground state wavefunction in first and second 
quantization: 

In first quantization, the ground state trial wavefunc- 
tion of BLQH in d — > limit is Halperin's (111) stated: 



Wi N 2 



) - Zj)Y[(Wi - Wj) 
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(Al) 

where z and w are the coordinates in layer 1 and layer 2 
respectively. As explicitly written, there are N\ electrons 
in the top layer and -/V2 electrons in the upper layer. 

In the EPQFM picture 2 ^, the ground state EPQFM 
wavefunction was written in second quantization form: 



where M — N is the angular momentum quantum num- 
ber corresponding to the edge of the system. 

In the following, we will show that the EPQFM wave- 
function in the second quantization form of Eqn lA2l is 
equivalent to (111) wavefunction in Eqn lAll 

The wavefunction Eqn lA2l can also be interpreted as 
the pairing of an electron in one layer and a hole in an- 
other which leads to exciton condensation. It is easy 
to show that the charges on the top layer are iVi = 
N cos 2 6/2 — Nui, while charges on the bottom layer 
are N 2 = N sin 2 6/2 = Nv 2 . Obviously Ni + N 2 = N. 

We can expand Eqn lA2l into: 
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and then transforming into the first quantization form 
leads to: 
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where A stands for anti-symmetrization. 

Moving the summation over the orbital states into A, 
the sum is essentially the (111) state: 
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This concludes that the orbital part of the projection of 
the EPQFM wavefunction EqnEHonto (Ni,N 2 ) sector 
is exactly the same as (111) state. If putting 6 — tt/2, we 
recover the balanced case in-2,. It is clear that the (111) 
wavefunctions with different (Ni, N 2 ) provide a complete 
basis, the ground state wavefunction Eqn lA2l can be ex- 
panded in terms of this basis, the angle 6 controls the 
expansion coefficients. 6 — tt/2(^ tt/2) corresponds to 
the balanced (imbalanced ) case. 



2. The fractional charges of merons 

In the subscetion, we will calculate the fractional 
charges of the excitations above the the (111) state. 

The charge density on each of the layers in the ground 
state Eqn lA2l is: 
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where N = iVj + N 2 and P is the permutation of N 
variables. 

Projecting onto the state with N\ electrons in layer 1 
and N 2 electrons in layer 2 
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The wavefunction of a meron in the bottom layer is: 
MS) 
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The charge density on each layer is: 
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FIG. 11: The charge densities of the smallest meron wave- 
function 



Subtracting Eqn lA8l from Eqn IAlOl leads to the frac- 
tional charge of the meron: 

Sp m i(z) = 0, Sp m2 (z) = -ia>[|^oO)| 2 ~ \<Pm{z)\ 2 ] 

( AU) 

The charges in top layer remains untouched, while v 2 
charge is moved from the origin to the boundary ( Fig.l 
) , namely, there is a hole with charge v 2 at the origin and 
extra charge —v 2 at the boundary. 

As concluded ir>22, because the meron wavefunction ig- 
nores the strong interlayer correlation, the charge distri- 
butions will be strongly modified by the neutral gapless 
mode. However, as shown in III B (2), the total charge 
is a topological property. Although the meron wavefunc- 
tion may not even be a qualitatively a good one, it still 
should give a correct answer on the total charge. Indeed, 
Eqn. IA1 II gives the correct total charge v 2 listed in Table 
2. 



APPENDIX B: QUANTUM HALL STATE TO 
WIGNER CRYSTAL TRANSITION IN SINGLE 
LAYER QH SYSTEMS 

Despite the success of wavefunction approach in SLQH, 
the nature of quantum phase transitions from QH to in- 
sulating state is beyond the scope of the wavefunction 
approach. It remains one of the outstanding problems 
in QH. Field theory is a must to study the transitions. 
Quantum Hall transitions in periodic potential in the 
presence of Coulomb interaction was investigated inSi in 
terms of Dirac fermion. The effects of both Coulomb in- 
teraction and disorders on the transition were discussed 
ir>Z!. But these investigations may not be directly rele- 
vant to the real QH transitions which happen in a con- 
tinuous system. In this appendix, inspired by the ESF 
to PSDW transition investigated in Sec. IV, we study the 
QH to Wigner Crystal (WC) transition in SLQH. 

It is instructive to point out the difference between 
SLQH and BLQH. In the SLQH, there is no true sym- 
metry breaking, the collective mode at k = turns out to 
be a local maximum. There is still a magneto-roton min- 
imum at qo ~ I (Fig.Bl), the collapsing of the minimum 
also signifies the collapsing of the QH gap, the system 
gets to a possible Wigner crystal state (Fig.Bl). Due to 
lack of true symmetry breaking, it is not known how to 
characterize the QH order by a local order parameter. 
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FIG. 12: The zero temperature phase diagram in a single 
layer quantum Hall system as the filling factor v changes. 
QH stands for Quantum Hall state, WC stands for Wigner 
Crystal state. The dispersion relations of the collective modes 
in the two phases are also shown. The cross in the WC means 
the negative magneto-roton minimum is replaced by the WC 
formation. 



The QH states may possess possible topological orders 
which need to be characterized by non-local order param- 
eters. So QH and WC also have two completely different 
orders: topological order and translational order. In the 
BLQH, the pseudo-spin sector is a charge neutral sec- 
tor, as shown ir*2&, both low energy modes at k = and 
k = ko are important. The collapsing of the magneto- 
roton minimum in the charge neutral sector leads to the 
PSDW in Fig.5. 

At fixed electron density, as the magnetic field in- 
creases, the filling factor v decreases, strong Quantum 
Hall states are discovered in v = 1/3,1/5,1/7, but not 
in 1/9 where the system may be in the WC state. The 
WC state will be pinned by disorders and become an in- 
sulating state. The quantunm Ginsburg-Landau action 
to describe the QH to the WC transition is: 

C wc {8p] = -Sp(-q 7 -u)[A p u 2 + r + c(q 2 - q 2 ) 2 ]5p(q,u) 



w{Sp) 3 +u(Sp) 4 



(Bl) 



where p — c^c is the electron density, Sp = p — po is the 
deviation from the average electron density, the momen- 
tum and frequency conservations in the cubic and quartic 
terms are assumed. 

In sharp contrast to the QGL action to describe the 
ESF to the PSDW transition Eqn [5^1 where the possible 
cubic interaction term is forbidden by the Z 2 exchange 
symmetry between the two layers, the cubic term is al- 
lowed. So Eqn IBll is similar to that EqnJSl] describing 
the SF to NS transition in Helium 4, except the former 
is in 2 + 1 dimension and the latter is in 3 + 1 dimen- 
sion. Unfortunately, the QH order is not transparent in 
Eqn|BlJ 

In Eqn IBll r is the magnetoroton minimum at q = qo 
which tunes the 1st order transition from the QH to the 
WC. In the QH fluid state, r > and p is uniform, < 
bp > = 0. In the WC, r < and < Sp >= £g n{G)e i3 - 3 
where / stands for the sum over all non-zero reciprical 
lattice vectors takes a lattice structure. Due to the cubic 
term, it can be shown that the most favorable lattice 
structure is the hexagonal lattice Fig. 6c. This is in sharp 
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contrast to the PSDW square lattice Fig. 5 in BLQH. 

Disorders will turn the 1st order transition into a 2nd 
order one, so scaling behaviors at low temperature are 
expected. As shown i n 76 ' 77 , Coulomb interactions will 
make the dynamic exponent of the 2nd order transition 
z = 1. We believe the Coulomb interaction is automat- 
ically incorporated in the action Eqn lBll How to incor- 
porate the effects of disorder into the action remains an 
outstanding problem. 



APPENDIX C: BILAYER QUANTUM HALL 
SYSTEM IN A PERIODIC POTENTIAL 

The Bilayer quantum Hall system in a periodic poten- 
tial can be described by the following relativistic model: 

£5 = \(du - ia,u)Zi\ 2 + m 2 \z\ 2 + u\zf + — l — e^xa^ax 

+ c/(M 2 -M 2 ) 2 (ci) 

wherei = 1, 2 are the two components, \z\ 2 = \zi\ 2 + \z2\ 2 , 
9 is the statistical angle and U > ( U < ) plays the 
role of the easy-plane ( Ising ) anisotropy which breaks 
the symmetry in the pseudo-spin sector from SU(2) down 
to U{1) x Z 2 . 

The SLQH case with just one complex field was inves- 
tigated ir>Z2. It can be considered as a model describing 
the transition from a quantum Hall state to a Mott in- 
sulator (MI) in a periodic potential as the strength of 
the periodic potential is varied. The fermionic version 
of the QH to MI transition in a periodic potential was 
investigated i n 76 i 77 . In both bosonic and fermionic ver- 
sions, the CS term is exactly marginal, the transition is a 
second order transition with continuously changing expo- 
nents depending on the statistical angle 9. One can view 
Eqn ICll to describe the transition in a BLQH system in 
a two 2-dimensional periodic potentials as the distance 
between the two layers is varied. We assume there is a 
critical distance d c such that m 2 (d < d c ) < 0,m 2 (d = 
d c ) = 0, m 2 (d > d c ) > 0. Due to the exact marginality of 
the CS term, we expect the transition is still 2nd order 
with continuously changing exponents depending on the 
statistical angle 9. 

In the imbalanced case, we add a Zeeman-likc term to 
EqniCT] 

C im b = A - MN 2 - N 2 ) (C2) 

The advantage of the relativistic theory Eqn lC2l over 
the non-relativistic model Eqn l23l in the main text is that 
it puts the fluctuations of spin sector and charge sector 
at the same footing. Of course, its physics could be dif- 
ferent from the real BLQH studied in the main text, but 
it is also interesting in its own right. It would be also in- 
teresting to compare its phases and phase diagram with 
those of the real BLQH discussed in the main text. In 
the following, we will first work out the phase diagrams 
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FIG. 13: The phase diagram of Eqn[C3] at T = (a) Easy- 
plane limit U > is self-dual. The three phases: XY Ferro- 
magnet ( F — XY ), exotic paramagnet ( P* ) and param- 
agnet ( P ) meet at the Tetra- quantum critical point which 
is a deconfined quantum critical poinlj 7 -^. The P* phase has 
a gapless photon mode, (b) Ising limit U < 0, the transition 
from the Ising ordered phase to the P* may be first order 
driven by gauge field fluctuations indicated by a dot. 

of two closely related models with non-compact and com- 
pact Maxwell gauge theories in (a) and (b), then we study 
the phase diagram of Eqn lC2l in (c). 

(a) Two component non-compact Maxwell theory 
If replacing the CS term in Eqn lC2l bv a non-compact 
Maxwell term, we get 

£non = \(d^-ia^)z l \ 2 +m 2 \z\ 2 + u\z\ i + ^(f^) 2 
+ U{\ Zl \ 2 -\z 2 \ 2 ) 2 -h z {\ Zl \ 2 -\z 2 \ 2 ) (C3) 

where / M „ = d^a^ - d v a^. 

This model had been studied in 78 in both isotropic 
U = and easy plane U > limit. Here we also study 
the Ising limit U < 0. In the easy-plane limit U > 0, 
the model was shown to be self-dual and there is a 2nd 
order transition through a deconfined quantum critical 
poirrfci 8 -. Under the duality transformation, the spin and 
charge sectors exchange their roles on the two sides of 
the QCP. There is one QPT driven by the imbalance on 
both sides of the QCP which are XY and inverted XY 
transitions respectively, so the two transitions on the two 
sides are indeed dual to each other ( Fig.Cla ) 78 . In the 
Ising limit U < 0, on the P* side, the transition is still in 
the universality class of the inverted XY transition, but 
in the Ising order side, there is no transition, the theory 
is not self-dual anymore ( Fig.Clb ). The transition may 
even be a 1st order transition driven by the gauge field 
fluctuations. 

The finite temperature phase diagram at the balanced 
case h z = is shown in the Fig. 2. There is one QPT 
driven by the temperature on both sides of the QCP 
which is in the universality class of the KT transition 
( Fig.C2a In the Ising limit U < 0, on the P* side, 
the transition is still in the universality class of the KT 
transition, but in the Ising ordered side, the transition is 
in the universality class of the Ising transition, the theory 
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FIG. 14: The finite temperature phase diagram of of Eqn lC3l 
at h z = 0. (a) Easy-plane limit U > which is self dual, (b) 
Ising limit U < which is not self dual 



FIG. 16: The finite temperature phase diagram of Eqn lC4l 
(a) Easy-plane limit U > (b) Ising limit U < 0. 
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FIG. 15: The phase diagram of the EqnE! at T = 0. (a) 
Easy-plane limit U > (b) Ising limit U < 0. 

is not self-dual anymore. ( Fig.C2b ). 

(b) Two component Compact Maxwell theory 
However, if the Maxwell term in Eqn lC3l is compact, 
then the compact gauge field confine the z a fields into a 
single unit vector n = z^(a) a f}Zf}, it can be shown that 
the compact model becomes: 

C comp = ^-{d^n) 2 + U(n z ) 2 - h z n z (C4) 

In the easy plane limit, the model is dual to a one 
component superconductor in an external magnetic field 
h z as described by Eqn[63l So there is one ( two ) QPT(s) 
for type I ( type II ) superconductor on the QD side only 
driven by the imbalance. The type II case was already 
shown in Fig. 9. The type I case was shown in Fig.C3a. 
The zero density transition was discussed in^i, it has the 
exact exponent z = 2,v = l,rj = 0. In the Ising limit, 
there is no transition on both sides as shown in Fig.C3b. 

The finite temperature phase diagram at the balanced 
case h z = is shown in the Fig.C4. On the ordered side, 
there is a KT transition in the easy-plane limit ( Fig.C4a 
) and an Ising transition in the Ising limit ( Fig.C4b ). 

So the phase diagram of the two component Maxwell 



theory is very different depending on the Maxwell term 
is compact or non-compact. 

(c) Two components Chern-Simon theory 
In the present model Eqn lC2l with the CS term, it was 
shown in 80 that in a pure compact Maxwell-Chern-Simon 
theory, monopolcs are linearly confined by a string of 
magnetic flux due to the CS term. Following the proce- 
dures inS, it is easy to see that in the model Eqn lC2l 
where the two component relativistic bosonic matter 
fields z a ,a — 1,2 are coupled to a compact CS term, 
the monopoles remain confined, therefore the two com- 
ponent bosonic fields are deconfined. So the CS term 
is compact or not in Eqn lC2l will not lead to qualita- 
tive changes in the phase diagram. In the balanced case 
h z = 0, in the easy-plane limit, due to the Zi symme- 
try, Z\ and Z2 have to condense at the same time. When 
d < d c ( m 2 < ), the system is in the ILCQH which is a 
Quantum Hall state with v = 1 integer Hall conductance 
o~xy = e 2 /h with a charge gap and a superfluid state with 
gapless Goldstone mode in the spin sector. While d > d c 
( m 2 > ), the system is in a quantum disordered state 
with a spin gap, but charge gapless. In the imbalanced 
case h z > 0, because z% and zi are deconfined, they do 
not have to condense at the same distance anymore. In 
the following, we discuss the Easy plane and Ising limit 
respectively. (1) Easy-palne limit, U > ( Fig.C5a ): On 
the QD side, when h z > to 2 , z\ becomes condensed, the 
QPT is in the same universality class as the QH transition 
in a single layer system, namely, one component complex 
scalar coupled to a CS gauge field in 2 + 1 dimension. 
Let's call this transition: QD-SLQH. Starting from the 
SLQH at fixed h z , if decreasing the distance further, the 
Z2 will also condense, the system will get into the BLQH 
state with interlayer phase coherence, let's call this tran- 
sition: SLQH-BLQH, because the U(l) CS gauge field 
is already broken by the z\ condensation, this transition 
is in the 3d XY universality class as shown in Fig.C5a. 
We can look at the transition from a different way: on 
the ILCQH side, when h z > — m 2 , z^, becomes uncon- 
densed. Obviously, the two transitions are not self-dual, 
in contrast to the non-compact Maxwell case in Fig.l. 
The QCP at d = d c in Fig.C5a is essentially a Tetra- 
quantum-critical point. (2) Ising limit U < ( Fig.C5b 
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). At the balanced case, on the ordered side, there is only 
a SLQH state which has no gapless Goldstone mode, so 
there is no BLQH state. At finite h z , there is a QD- 
SLQH transition in the disordered side, while there is no 
transition in the SLQH side. Note that the transitions in 
both easy plane limit (a) and Ising limit (b) are second 
order due to the exact marginality of the CS term, even 
the theory is not self-dual in both cases. This is in sharp 
contrast to the corresponding phase diagram Fig. CI in 
the non-compact Maxwell theory. 



FIG. 17: The phase diagram of of Eqn[C2]at T = 0. (a) Easy- 
plane limit U > 0, there are three phases meeting at the Tetra- 
quantum critical point: (1) Bilayer Quantum Hall state with 
interlayer coherence < Z\ >^ 0, < ( BLQH+ILCH ) 

(2) Single Quantum Hall state < Z\ 0, < z 2 >= (SLQH 
) (3) Quantum Disordered state < Zj > — < Z2 >= 0. (b) 
Ising limit U < 0, there is no BLQH state. The SLQH and 
the QD meets at a Bi-quantum critical point. 
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FIG. 18: The finite temperature phase diagram of Eqn lC2 
(a) Easy-plane limit U > (b) Ising limit U < 



The finite temperature phase diagram at the balanced 
case h z = is shown in the Fig.C6. On the ordered side, 
there is a KT transition in the easy-plane limit ( Fig.C6a 
) and an Ising transition in the Ising limit ( Fig.C6b ). 
Fig.C6 is similar to Fig.C4 to some extent. 



One can compare Fig.C5(a) ( Fif.C6(a) ) with Fig. 9(b) 
( Fig. 8 ) in the main text. As expected, they arc different, 
because they are two different systems. The former is the 
BLQH system in a periodic potential, while the latter is 
the real BLQH in a continuous system. 
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